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Abstract 

This  report  presents  the  results  of  initial  testing  of  an  actively  blown  grid  for  producing 
high-intensity,  large  length-scale  freestream  turbulence  for  use  in  a  linear,  transonic 
turbine  cascade.  The  design  was  achieved  using  data  presented  in  the  work  of  several 
other  researchers.  The  desired  turbulence  intensity  (Tu%)  and  integral  length  scale 

(Ax)  to  match  gas  turbine  engine  conditions  with  similarity  are  Tu  =  15-20%  and  A*  ~  4 
cm.  The  testing  stage  of  this  work  was  performed  in  a  low  speed  test  section  of  similar 
flow  area  as  the  turbine  cascade.  Data  are  presented  for  the  grid  both  with  and  without 
active  blowing.  The  grid  data  without  blowing  was  compared  to  correlations  of 
turbulence  and  length  scale  decay  presented  by  Roach^  and  Baines  and  Peterson’  with 
good  agreement  for  turbulence  intensity  decay  and  fair  for  length  scale  dispersion.  Initial 
tests  achieved  a  turbulence  intensity  of  approximately  17%  and  a  length  scale  of 
approximately  3.5  cm  at  a  distance  of  18  bar  widths  downstream  of  the  blown  grid. 


Introduction 

Using  design  ideas  from  published  Uterature^  ’’^  '’’^,  an  actively  blown  turbulence 
generating  grid  was  developed  incorporating  both  upstream  and  downstream  blowing  to 
generate  large  length-scale,  high  intensity  freestream  turbulence  for  use  in  a  6”  span 
transonic  turbine  cascade.  The  generator  (shown  in  Figure  1)  has  been  designed  to 
produce  freestream  turbulence  with  a  turbulence  intensity  (Tu%)  of  approximately  15- 
20%  with  an  integral  length  scale  of  4  cm.  This  turbulence  level  (intensity  and  scale)  was 
determined  through  research  done  at  NASA  Glenn  by  Dr.  Jim  Van  Fossen,  et  al®.  Van 
Fossen  used  a  60°  section  of  a  GE  90  combustor  in  a  wind  tunnel  to  model  the  exit  flow 
of  a  combustor,  including  inlet  swirl  vanes,  film  cooling  holes  and  dilution  holes.  AU 
three  of  these  combustor  components  contribute  strongly  to  the  combustor  exit  turbulence 
conditions.  Hot  wire  surveys  were  performed  downstream  of  the  combustor  to 
characterize  the  turbulence  conditions  at  the  exit.  Combustion  was  not  included,  but 
previous  research  by  Moss  and  Oldfield^  showed  that  turbulence  intensity  and  length 
scale  were  not  significantly  affected  by  the  presence  of  combustion.  Using  Van  Fossen’s 
results,  the  desired  length  scale  was  increased  to  match  the  chord  length  of  the  blades 
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Note:  All  units  are  in  inches. 


Figure  1:  Schematic  of  Actively-Blown  Turbulence  Generating  Grid 


used  in  the  Virginia  Tech  transonic  turbine  cascade  facility  with  similarity  to  the  blades 
of  a  GE  90  engine.  Van  Fossen  reported  length-scales  between  1  and  1.5  cm  from  his 
combustor  simulation,  where  the  chord  length  of  a  first  stage  GE  90  turbine  blade  is 
estimated  as  approximately  2  in  (5.08  cm).  Our  turbine  blade  has  a  chord  length  of  5.4  in 
(13.72  cm),  so  the  turbulent  length  scale  desired  was  determined  to  be  approximately  4 
cm.  Testing  of  the  grid  was  performed  in  a  low  speed  environment  to  characterize  the 
turbulence  from  the  grid.  The  mass  flow  ejected  from  the  grid  holes  was  adjusted  in  an 
attempt  to  reach  the  turbulence  intensity  and  length  scale  desired.  The  momentum  flux 
ratios  (blowing  vs.  freestream)  can  then  be  matched  in  the  transonic  facility  to  hopefully 
achieve  the  same  results  in  both  environments  (low-speed  versus  transonic). 

The  test  setup  for  the  low-speed  tests  is  shown  in  Figure  2.  The  setup  includes  a  5  HP 
blower  attached  to  several  transition  pieces,  followed  by  a  test  section.  The  blower 
created  mainstream  flow  velocities  of  approximately  20  m/s,  whereas  the  inlet  velocity  of 
the  transonic  wind  tunnel  is  approximately  100  m/s.  The  first  transition  piece  changes 
the  cross  section  of  the  flow  path  from  circular  to  rectangular,  the  second  transition  piece 
changes  the  area  to  accommodate  the  turbulence  generating  grid  and  the  third  transition 
changes  the  flow  area  from  the  cross  section  of  the  grid  to  the  cross  section  of  the 
transonic  wind  tunnel  test  section,  which  is  12  inches  high  by  6  inches  wide.  All 
transition  pieces  were  designed  to  avoid  flow  separation.  A  Plexiglas  test  section,  which 
has  this  same  cross-sectional  geometry  and  is  20  inches  in  length,  was  fabricated  to  take 
data  downstream  of  the  turbulence  grid  in  the  low-speed  environment.  Initial  testing  of 
the  blower  through  downstream  hotwire  velocity  traverses  indicated  that  the  background 
freestream  turbulence  levels  from  the  blower  were  very  high  (approximately  12%).  To 
lower  the  background  level  of  turbulence  entering  the  test  section,  a  flow  straightener 
section  was  installed  which  consists  of  a  section  of  honey  comb,  followed  by  a  fine  mesh 
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Figure  2:  Test  Setup  for  Low  Speed  Testing  of  Turbulence  Grid 


screen,  then  a  straight  section  12  inches  in  length  followed  by  another  section  of 
honeycomb.  This  flow  straightener  was  installed  just  upstream  of  the  rectangular 
transition  piece  in  Figure  2.  Background  turbulence  levels  were  decreased  to  slightly  less 
than  2%. 


Blowing  air  for  the  turbulence  generating  grid  is  supplied  through  high-pressing  lines  to  a 
plenum  tank.  The  plenum  then  supplies  air  to  each  of  the  three  blowing  bars  (each  square 
bar  is  1”  on  a  side)  of  the  grid.  A  diagram  of  the  blowing  supply  setup  is  shown  in  Figure 
3.  The  grid  has  1/8  in  diameter  holes  drilled  on  both  the  upstream  and  downstream  sides, 
spaced  1  inch  apart.  The  bar  spacing  was  set  to  3  inches  based  on  a  combination  of 
results  from  other  researchers  mentioned  earlier. 


Ball  Valves  to  Control 
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Figure  3:  EHagram  of  Blowing  Air  Supply  System 
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Velocity  measurements  (including  fluctuating  component)  were  taken  at  various 
locations  throughout  the  test  section  using  a  hot  wire  anemometer  probe.  The  probe  was 
manually  traversed  vertically  at  one  downstream  location  to  demonstrate  flow  uniformity 
(checking  if  wakes  are  present  from  bars  on  grid)  in  the  test  section.  This  traverse  was 
performed  at  13  bar  widths  downstream  of  the  grid.  Four  locations  in  the  streamwise 
direction  downstream  of  the  grid  were  chosen  to  take  data  in  the  center  of  the  passage  in 
the  absence  of  an  adequate  traversing  system  (to  be  added  soon).  Statistical  and  spectral 
analysis  of  the  velocity  data  at  these  four  locations  was  used  to  calculate  the  turbulence 
intensity  and  length-scale  and  the  decay/dispersion  of  these  quantities  as  the  flow 
progresses  through  the  t  jst  section.  The  data  were  sampled  at  20  kHz,  filtered  at  10  kHz 
and  AC  coupled  with  a  5V  offset  to  increase  dynamic  resolution.  All  data  points  were 
sampled  for  6  s  and  4  samples  were  taken  at  each  point  for  averaging  purposes. 

From  pressure  and  temperature  measurements  at  the  plenum  and  grid  blowing  hole  exit, 
which  have  not  yet  been  performed,  the  momentum  flux  ratio  (blowing/freestream)  can 
be  determined  for  the  low  speed  tests.  The  momentum  flux  ratio  will  be  useful  for 
matching  turbulence  characteristics  when  the  turbulence  grid  is  placed  in  the  flow  field  of 
the  transonic  tunnel  for  use  in  turbine  blade  beat  transfer  testing. 

Results 

Characterization  of  Turbulence :  Statistical  and  Spectral  Analysis  of  Velocity  Data 

In  order  to  characterize  the  turbulence  from  the  developed  actively  blown  grid,  statistical 
and  spectral  analyses  were  performed  on  all  data,  both  with  and  without  blowing 
activated.  Velocity  data  at  each  sampled  point  in  the  flowfield  was  actually  averaged 
over  6  data  samples  to  give  better  statistics  and  spectra.  Statistical  and  spectral  analyses 
included  calculation  of  the  mean  velocity  (U),  fluctuating  component  of  velocity  (u),  root 
mean  square  (RMS)  of  the  fluctuating  component  of  velocity  (m'),  autocorrelation  (Rv), 
probability  density  functions  (PDPO.  and  power  spectral  density  (PSD)  calculations. 

The  streamwise  instantaneous  velocity  signal  can  be  decomposed  into  a  mean  and 
fluctuating  component  a*^  follows; 


u  =U  +  u 


where  U  is  the  mean  velocity  and  u  is  the  fluctuating  component  of  velocity.  The  mean, 
fluctuating  component  and  RMS  of  the  discretized  velocity  data  were  calculated  as 
follows: 


U  = 


u=u  — 


=  u-U 


u 


xl/2 
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Using  these  values,  the  turbulence  intensity  (Tu%)  of  the  flow  at  a  given  point  was 
calculated  (assuming  isotropic  turbulence)  by: 


Tu%=  — 
U 


The  autocorrelation  (Ry)  of  the  fluctuating  component  of  velocity  (normalized  such  that 
Ry  is  equal  to  1  at  zero  lag,  T=0),  was  calculated  by: 


/?v(7)  = 


u{t)u{t  +  x) 

17^ 


N 


'•k  '' 


*+y 


i=l 


U 


n 


where  r  =  j  •  Af 


The  autocorrelation  was  used  to  determine  the  integral  time  scale  (3)  by  integrating 
under  the  curve  to  the  first  zero  crossing. 


3=  J/?,(T)dT=|i/?,.Ar 

f=0 

The  streamwise  integral  length  scale,  which  is  representative  of  the  largest  eddies  in  the 
turbulent  flowfield,  was  then  determined  by  invoking  Taylor’s  hypothesis  of  frozen 
turbulence: 


Ax  =  U.3 

where  No  is  the  point  of  the  first  zero  crossing. 

The  probability  density  function  (PDF)  of  each  set,  of  velocity  data  was  also  determined. 
The  PDF  is  useful  in  visualizing  the  compar.jOn  of  mean  velocity  and  deviation  of  the 
data  from  the  mean  (proportional  to  Tu%)  of  different  data  sets.  The  PDF  also 
demonstrates  the  Gaussian  behavior  of  the  data.  The  PDF  of  a  data  set  was  done  by 
“binning”  the  velocity  data  into  discrete  bins  (histogram)  and  normalizing  by: 

N- 

PDF  = - - 

^Tou^  -  Am 


Finally,  the  power  spectral  density  (PSD)  each  data  set  was  determined  by: 

^^^_2FFT(u)FFT\u)  .  . . 

- — -  where  *  denotes  the  complex  conjugate 


The  power  spectral  density  was  frequency  averaged  to  give  a  cleaner  curve  by  averaging 
PSD  data  over  a  fi:equency  range  of  Af=l(X)  Hz  (data  points  over  every  100  Hz  were 
averaged  to  be  a  single  point).  The  power  spectral  density  gives  an  indication  of 
turbulent  energy  content  of  the  flow  at  different  frequencies,  which  are  inversely 
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proportional  to  the  size  of  the  turbulent  eddies  in  the  flow  (the  wave  number,  K,  is  equal 
to  2t[S\J).  For  comparison,  in  the  discussion  section,  the  PSD  data  was  normalized  by: 


U 


and  compared  to  the  Von  Karman  spectrum  equation: 


E(f)U 

A, 


=  4- 


v-5/6 


Using  these  statistical  and  spectral  analysis  tools,  the  characteristics  of  the  turbulence 
could  be  determined. 


Experimental  Turbulence  Data 

This  section  uses  the  tools  described  to  characterize  the  turbulence  from  the  designed  grid 
experiments  both  with  and  without  active  blowing.  As  discussed,  four  locations  spaced 
in  the  streamwise  direction  were  chosen  to  take  hotwire  measurements  at  mid-passage  (3 
inches  from  the  sidewalls  and  6  inches  from  the  top  and  bottom  walls).  Mid  passage  is 
directly  behind  one  of  the  blowing  bars  of  the  grid.  Table  1  gives  the  axial  location  in  bar 
widths  (bars  are  1  inch  in  width)  downstream  of  the  grid: 


9.5 

13.5 

17.5 

21.5 

Table  1:  Locations  of  Data  Acquisition  Points 


The  hotwire  data  at  each  of  these  points  was  reduced  using  a  Matlab  code  to  calculate 
statistics  and  autocorrelation  to  determine  turbulence  intensity  and  integral  length  scale. 
A  sample  of  the  Matlab  code  file  written  to  accomplish  all  calculations  in  this  report  is 
included  as  Appendix  A.  The  mean  flow  velocity,  turbulence  intensity,  and  integral 
length  scale  at  each  location  for  different  blowing  rates  from  the  bars  (including  no 
blowing)  are  shown  in  Figure  4  through  Figure  6,  respectively.  In  the  absence  of  a 
system  of  measuring  plenum  pressure  independently  for  each  bar,  the  blowing  rate  was 
estimated  by  how  far  the  upstream  ball  valves  between  the  blowing  air  plenum  and  grid 
bars  were  opened.  There  is  some  run-to-run  error  inherent  in  this  method. 
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Plot  of  Mean  Velocity  at  4  Streamwise  Positions  Behind  Center 
Bar  at  Different  Blowing  Rates 
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Figure  4.  Plot  of  Mean  Velocity  at  4  Streamwise  Locations  at  .Different  Blowing  Rates 


Plot  ofTurbulence  Intensity  at  4  Streamwise  Positions  Behind 
Center  Bar  at  Different  Blowing  Rates 
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Figure  5.  Plot  of  Turbulence  Intensity  at  4  Streamwise  Locations  at  Different  Blowing  Rates 
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Plot  of  Integral  Length  Scale  at  4  Streamwise  Positions  Behind  Center 
Bar  at  Different  Blowing  Rates 
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Figure  6.  Plot  of  Integral  Length  Scale  at  4  Streamwise  Locations  at  Diffa-ent  Blowing  Rates 

The  mean  flow  in  the  test  section  without  the  grid  installed  was  approximately  20  m/s  at 
center  passage,  which  is  independent  of  streamwise  measurement  location.  Figure  4 
demonstrates  that  at  locations  near  the  grid  without  blowing,  there  is  a  velocity  deficit 
due  to  the  wake  of  the  grid  bars.  This  velocity  deficit  decreases  as  you  march  in  the 
streamwise  direction  and  away  fi-om  the  influence  of  the  wake.  The  general  trend  in 
velocity  with  blowing  is  that  the  velocity  at  each  location  increases  with  increased 
blowing  as  the  bar  wake  is  filled  in  with  higher  momentum  fluid  and  eventuaUy  a  jet 
forms.  The  initial  decrease  in  mean  velocity  at  the  lowest  blowing  rate  at  Location  1 
could  not  be  explained.  The  decrease  in  mean  velocity  at  the  highest  blowing  rate  is  due 
to  the  plenum  pressure  dropping,  most  likely  due  to  the  inlet  flow  to  the  plenum  choking, 
rhe  turbulence  intensity  of  the  flow  without  the  grid  installed  was  approximately  2%. 
Figure  5  demonstrates  that  with  the  grid  in  place  without  blowing,  the  intensity  increases 
significantly  to  approximately  10-15%  depending  on  streamwise  location.  As  blowing  is 
activated  and  increased,  there  is  first  a  decrease  in  turbulence  intensity,  then  a  steady 
increase  with  blowing  rate.  The  initial  decrease  is  probably  due  to  filling  in  the  wake  of 
the  bars,  which  essentially  reduces  the  effectiveness  (or  effective  bar  width)  of  the  grid. 
The  integral  length  scale  of  the  flow  was  shown,  with  the  exception  of  one  point,  to 
steadily  increase  with  blowing  rate  as  can  be  seen  in  Figure  6.  The  length  scale  at  higher 
blowing  rates  is  seen  to  “settle  out”  with  no  further  increase.  For  the  data  presented  here, 
the  blowing  rate  (blowing  plenum  pressure)  is  Just  an  initial  setting.  The  current  results 
show  that  at  approximately  18  bar  widths  downstream  of  the  grid,  the  turbulence  intensity 
is  as  high  as  17%  with  an  integral  length  scale  of  3.5  cm.  Continuing  work  is  being  done 
by  myself  to  determine  the  optimal  plenum  pressure  and  internal  bar  pressure  that  will 
provide  the  desired  turbulence  intensity  and  length  scale  at  18  bar  widths  downstream  of 
the  grid,  as  well  as  investigate  the  uniformity  of  the  flow  from  the  blown  grid  in  the 
cross-stream  direction. 
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Statistics  and  Power  Spectra 

At  each  measurement  location,  the  PDF,  autocorrelation  and  PSD  were  determined  and 
the  results  compared  for  each  blowing  rate.  Figure  7  through  Figure  9  present  the  PDFs, 
autocorrelations,  and  PSDs  at  each  location  for  varying  blowing  rates,  respectively. 


Figure  7.  Probability  Density  Functions  (PDF)  at  Each  Streamwise  Measurement  Location  at 

Different  Blowing  Rates 


The  probability  density  functions  at  different  blowing  rates  for  each  measurement 
location  (Figure  7)  demonstrate  that  as  blowing  is  increased,  the  mean  velocity  and 
deviation  from  the  mean  (which  is  proportional  to  the  turbulence  intensity)  change  in  the 
manner  discussed  in  the  previous  section.  With  the  exception  of  one  point  at  Location  I, 
the  mean  velocity  tends  to  increase  and  the  turbulence  intensity  (“thickness”  of  the  PDF) 
tends  to  increase  as  the  blowing  rate  is  increased.  The  PDFs  demonstrate  the  Gaussian 
behavior  of  the  velocity  data. 

Through  observation  of  the  autocorrelation  curves  (Figure  8),  it  can  be  seen  that  the  first 
zero  crossing  value  does  not  change  significantly  with  increased  blowing  at  each 
streamwise  location.  However,  the  shape  of  the  curves  changes  with  increased  blowing, 
so  integrating  under  the  curve  to  obtain  the  integral  time  scale  will  yield  increasing  time 
scales  with  increasing  blowing  rate.  This  fact,  coupled  with  the  mean  velocity  change 
will,  in  general,  cause  a  corresponding  increase  in  length  scale  with  increasing  blowing 
rate.  This  is  in  agreement  with  the  results  shown  in  Figure  6. 
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omarantBloHilng  RrtM 


Figure  8.  AutcxOTrelations  (Ry)  at  Each  Stream  wise  Measurement  Location  at  Different  Blowing  Rates 


Figure  9.  Power  Spectral  Density  (PSD)  at  Each  Streamwise  Measurement  Location  at  Different  Blowing 

Rates 
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The  power  spectral  density  (PSD)  plots  for  each  location  confirm  the  results  for  the 
turbulence  intensity.  As  mentioned  earlier,  power  spectral  density  gives  an  indication  of 
turbulent  energy  content  of  the  flow  at  different  frequencies,  which  is  inversely 
proportional  to  the  size  of  the  eddies.  The  PSD  plots  for  each  location  at  different 
blowing  rates  (Figure  9)  demonstrate  that  as  blowing  is  increased  the  turbulent  energy  at 
each  frequency  (or  wave  number)  increases.  Integrating  under  the  PSD  will  yield  the 
RMS  of  the  velocity  fluctuations  («'),  which  is  proportional  to  turbulence  intensity.  At 
each  location,  the  PSD  without  blowing  can  be  seen  to  have  a  “hump”  around  100  Hz. 
This  is  believed  to  be  a  product  of  weak  vortex  shedding  from  the  bars  without  blowing. 
This  phenomenon,  and  a  more  detailed  analysis  of  the  PSD  for  sel^^t  data  sets;  will  be 
discussed  in  detail  in  the  discussion  section. 

Discussion 

Comparison  to  Grid  Generated  Turbulence 

The  set  of  tests  performed  involved  collecting  data  at  four  streamwise  positions  to 
compare  the  published  data  to  grid  generated  turbulence  studies.  The  “grid”  in  this  work 
is  actually  a  series  of  horizontal  bars  rather  than  a  2-D  grid  arrangement,  however,  it  is 
suggested  that  the  data  should  stiU  match  well  with  grid  generated  turbulence 
correlations.  Studies  performed  by  Roach^  and  Baines  and  Peterson^  demonstrated  the 
streamwise  variation  of  turbulence  intensity  and  integral  length  scale  behind  tube  bundles 
and  screens.  From  these  studies,  it  has  been  shown  that  the  turbulence  decay  and  length 
scale  dispersion  can  be  characterized  by  equations  in  the  form; 


Tu%  =  c- 


and  K=ad- 


where  c  and  a  are  constants  given  as  0.80  and  0.20,  respectively,  x  is  the  streamwise 
distance,  d  is  the  bar  width,  and  n  and  m  are  exponents.  The  value  of  n  is  typically  -5/7 
and  m  is  between  0.5  and  0.56  according  to  the  published  data.  The  daia  collected  behind 
the  grid  in  this  work  was  compared  to  the  published  curves.  Figure  IP  through  Figure  12 
show  the  mean  velocity,  turbulence  intensity,  and  integral  length  scale  versus  streamwise 
measurement  location  both  with  and  without  blowing.  One  blowing  rate  was  chosen  for 
simplicity,  the  “valve  half  open”  condition.  From  Figure  10,  it  can  be  seen  that  without 
blowing,  the  velocity  data  increases  slightly  as  you  march  downstream  and  get  out  of  the 
influence  of  the  bar  wake,  as  discussed  earlier.  The  same  type  of  trend  is  seen  for  the 
velocity  with  blowing,  however,  the  velocity  increases  much  faster  as  a  result  of  the 
blowing  filling  in  the  wake.  The  streamwise  decay  of  turbulence  intensity  is  shown  in 
Figure  1 1.  The  data  with  and  without  blowing  is  plotted  along  with  the  turbulence  decay 
equation  given  by  Roach  and  Baines  and  Peterson.  Similarly,  the  streamwise  dispersion 
of  turbulent  length  scale  is  shown  in  Figure  12.  These  two  figures  demonstrate  that 
blowing  follows  a  different  “decay  law”  than  grid  generated  turbulence.  The  decay  of 
turbulence  intensity  without  blowing  agrees  well  with  the  correlations,  however  the 
integral  length  scale  does  not  agree  as  well  Deviation  from  the  correlations  could  be  due 
to  the  geometry  used  here  (Le.  horizontal  square  bars  versus  round  wire  mesh  grids).  It  is 
also  possible  that  the  discrepancy  is  due  to  a  contraction  being  present  downstream  of  the 
grid.  The  turbulence  intensity  data  with  blowing  was  shown  to  obey  a  power  law  of 
Tu%=1.5(x/d)°'^*\  however  this  was  only  valid  at  this  blowing  rate.  Increasing  or 
decreasing  blowing  rate  changed  this  relationship. 
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Plot  of  Mean  Velocity  vs  Streamwise  Location 
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Figure  10.  Variation  of  Mean  Velocity  in  the  Streamwise  Direction  with  and  without  Blowing 


Plot  of  Turtxilence  Intensity  vs  Streamwise 
Location  Compared  to  Grid  Generated  Turbulence 
Correlation 


Figure  1 1.  Variation  Turbulence  Intensity  in  the  Streamwise  Direction  with  and  without  Blowing 


Figure  12.  Variation  Length  Scale  in  the  Streamwise  Direction  with  and  without  Blowing 
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The  probability  density  fiinction  of  the  data  with  and  without  blowing,  shown  in  Figure 
13  demonstrate,  as  before,  the  increase  in  the  mean  velocity  and  turbulence  intensity  in 
the  streamwise  direction.  Again,  it  can  be  seen  from  the  PDFs  that  as  we  march  along  in 
the  streamwise  direction,  the  mean  velocity  increases  and  the  deviation  becomes  smaller, 
indicating  a  decrease  in  turbulence  intensity.  The  relative  widths  at  each  location  with 
and  without  blowing  indicate  that  the  turbulence  intensity  is  higher  for  the  blowing  case 
and  the  mean  increases  more  with  blowing  since  the  velocity  deficit  is  overcome  earlier 
due  to  the  jets  filling  in  the  wake  behind  the  bars.  Of  interest,  and  as  yet  unexplained,  is 
that  the  PDFs  for  the  unblown  cases  (and  in  the  blown  cases  to  a  lesser  extent)  collapse 
down  to  a  common  curve  on  the  high  side  (+3o)  of  the  mean  velocity. 


PDPk  «t  4  StwmwiM  Potitom  BMnd  Canlv  Bar  wtm  Btoning 


Figure  13.  Probability  Density  Functions  of  Velocity  fcff  Various  Streamwise  Positions  with  and 

without  Blowing 

Determination  of  Vortex  Shedding  for  Unblown  Grid 

In  the  results  presented  in  the  spectra  for  the  velocity  data,  it  was  indicated  that  the 
spectra  for  the  unblown  grid  data  displayed  a  “hump”  in  the  data  around  100  Hz.  The 
possibility  of  vortex  shedding  occurring  from  the  bars  was  investigated.  The  vortex 
shedding  frequency  of  a  body  can  be  described  by  the  dimensionless  parameter  of 
Strouhal  (S), 

/  D 
- 

U 

where  fs  is  the  vortex  shedding  frequency,  D  is  the  width  of  the  bar,  and  U  is  the  mean 
velocity  upstream  of  the  bar.  The  Strouhal  number  is  a  function  of  the  geometry  of  the 
body  and  Reynolds  number.  For  a  square  bar  in  cross-flow,  the  Strouhal  number  has 
been  demonstrated"’  "  to  be  0.13  (versus  0.21  for  circular  cylinders)  with  strongest  vortex 
shedding  in  the  range  of  Reo  ~  10^.  Using  a  value  of  S=0.13,  a  mean  velocity  of  19  m/s 
upstream  of  the  bar  and  a  diameter  of  1  inch  (0.0254  m),  the  shedding  frequency  was 
determined  to  be  97  Hz.  Figure  14  shows  the  power  spectral  density  at  each  streamwise 
location  without  blowing.  It  can  be  seen  from  the  plot  that  the  “hump”  in  the  PSD  is 
centered  on  approximately  100  Hz.  The  magnitude  decreases  with  streamwise  distance, 
suggesting  that  the  vortex  shedding  is  relatively  weak  and  the  vortices  break  up  as  the 
flow  progresses  downstream.  This  makes  sense,  since  the  Reynolds  number  based  on  the 
diameter  of  the  bar  is  only  ReD=32,000,  whereas  according  to  the  literature,  vortex 
shedding  for  square  bars  occurs  at  a  Rcd  ~  10^.  Therefore,  we  would  not  expect  strong 
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vortex  shedding  from  the  bars  in  this  velocity  regime.  However,  this  is  of  importance  in 
continuing  work,  as  vortex  shedding  from  the  bars  in  the  transonic  facility,  where  the 
velocity  at  the  bars  can  be  as  high  as  120  m/s  may  result  in  stronger  vortex  shedding. 
Blowing  should  eliminate  concerns  with  vortex  shedding,  but  this  analysis  is  worth 
consideration  in  the  transonic  tunnel. 


Frequency  Averaged  Power  Spectral  Density  at  4  Streamwise 
Positions  Behind  Center  Bar 


Frequency  (Hz) 


Figure  14.  Powa*  Spectra  at  Each  Measurement  Lx)cation  without  Blowing 


Comparison  of  Power  Spectral  Density  of  Turbulent  Flow  with  Von  Kdrmdn  Spectra 

For  purely  isotropic  turbulence  in  the  absence  of  outside  forces  and  viscous  dissipation, 
the  Von  Karman  spectrum  equation  holds: 


E{f)-U 

A, 


=  4- 


1  + 


■  A, 

3-U 


2  \ 


-5/6 


where  E(f)  is  the  power  spectral  density.  If  the  PSD  for  the  velocity  data  is  normalized 
by: 


and  plotted  versus 


2^ -A, 
U 


=  k-A^  where  k  is  the  wave  number 


the  results  can  be  directly  compared  to  the  Von  Karman  spectrum  equation. 

Figure  15  and  Figure  16  are  comparisons  of  two  PSDs  (one  with  blowing  and  one 
without  blowing  20  bar  widths  downstream  of  the  grid)  from  this  study  with  the  Von 
Karman  equation.  There  are  three  regions  of  note  in  the  PSD.  The  first  region  is  the 
relatively  flat  region  at  the  low  wave  numbers  at  the  beginning  of  the  curve,  which 
contains  the  relatively  low  energy,  large  length  scale  turbulent  eddies.  At  the  end  of  this 
region  is  where  the  most  energetic  eddies  exist,  which  occurs  where  k-Ax~1  (i.e.  the  wave 
number  is  approximately  equal  to  the  inverse  of  the  integral  length  scale).  This  wave 
number  is  denoted  as  K^.  Through  observation  of  a  plot  of  the  wave  number  and 
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normalized  power  spectra  product  (shown  in  Figure  17  and  Figure  18),  it  can  be  seen 
from  the  spectra  that  the  most  energetic  scales  should  occur  at  k-A*  ~1.  The  data  with 
blowing  match  up  fairly  well  with  the  Von  Karman  curves,  however,  an  interesting  point 
to  note  is  that  in  both  Figure  15  and  Figure  17  (data  for  grid  without  blowing),  there  is  a 
slight  departure  of  the  data  from  the  theoretical  Von  Karman  curves.  This  is  believed  to 
be  due  to  the  vortex  shedding  discussed  earher.  In  Figure  17,  using  the  previously 
calculated  vortex  shedding  frequency  of  97  Hz,  and  for  the  case  plotted  in  the  figures,  U 
was  equal  to  19.91  m/s  and  A*  was  equal  to  2.06  cm,  yields  k-Ax=  0.778.  This  line  is 
shown  in  both  plots.  This  is  close  to  k  A^  ~1,  but  as  can  be  seen,  there  is  a  much  higher 
peak  in  the  PSD  of  the  unblown  grid  comps^ied  to  the  Von  Karman  curve.  It  is  possible 
that  coherent  vortex  shedding  structures  are  still  in  existence  20  bar  widths  downstream 
of  the  grid,  although  they  do  not  seem  to  be  very  strong. 

The  next  region  is  a  region  where  the  spectrum  foUows  a  power  law  decay.  This 
region  is  also  caUed  the  inertial  subrange  and  is  the  region  where  the  higher  frequency 
(smaller  scale)  turbulence  exists  independently  of  larger  scales  and  energy  is  cascaded 
down  from  large  scales  to  small  scales  until  the  eddies  are  sufEiciently  small  enough  to  be 
m  the  range  of  viscous  dissipation.  This  region  is  useful  in  determination  of  the 
dissipation  rate,  e,  of  the  turbulence.  By  curve-fitting  the  data  in  the  inertial  subrange  to 
an  equation  of  the  form  E(K)=A  e^^  k'^^,  the  dissipation  rate  can  be  determined.  The 
value  of  A  varies  depending  on  the  source,  with  values  of  0.50-0.53  being  reported 
The  dissipation  rate  was  not  calculated  for  this  work.  The  third  region  is  the  departure 
from  the  -5/3  decay  law  and  is  a  product  of  viscous  dissipation  and  occurs  at  very  small 
scales  (Kolmogorav  or  dissipation  length  scales).  This  region  is  readily  evident  in  both 
the  unblown  and  blown  grid  data.  The  inertial  subrange  of  the  blown  grid  seems  to  be 
much  larger  than  that  of  the  unblown  grid,  as  the  data  for  the  blown  grid  holds  to  the  Von 
Karman  curve  and  -5/3  slope  for  a  larger  range  of  wave  numbers. 


Frequency  Averaged  Power  Spectral  Density  without  Blowing  Compared 
to  Von  Karman  Cpectrum  Equation 


Figure  15.  PSD  Data  without  Blowing  Compared  to  Von  Karmen  Relation 
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Frequency  Averaged  Power  Spectral  Density  with  Blowing  Compared  to 

Von  Karman  Spectrum  Equation 
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Figure  16.  PSD  Data  with  Blowing  Conqjared  to  Von  Kdrmdn  Relation 


Energy  Content  at  Given  Wave  Number  without  Blowing  Compared  to 
Von  Karman  Spectrum  Equation 


Figure  17.  Energy  Content  versus  Wave  NumbCT  without  Blowing  Compared  to  Von  K^mdn 

Relation 
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Energy  Content  at  Given  Wave  Number  with  Blowing  Compared  to  Von 
Karman  Spectrum  Equatton 


2*K‘AJU  or  ic*A, 


Figure  18.  Energy  Content  versus  Wave  NumbCT  with  Blowing  Compared  to  Von  K^miin 

Relation 

The  wave  number,  Kj,  where  viscous  dissipation  starts  can  be  related  to  the  Kolmogorav 
length  scale  (Tj)  by: 


0.1 


The  value  of  can  be  estimated  by  plotting  the  dissipation  energy  spectrum  K^-  E(f). 
Figure  19  shows  the  dissipation  spectrum  both  with  and  without  blowing.  The  value  of 
Kd  is  defined  as  the  maximum  of  these  curves,  where  the  derivative  of  the  dissipation 
spectrum  with  respect  to  wave  number  is  zero.  From  comparison  of  the  curves,  it  is  not 
readily  apparent  which  has  a  higher  Kd.  However,  these  curves  do  provide  some  insight 
as  to  the  relative  size  of  the  Kolmogorav  length  scale,  within  the  limitations  of  the 
measurement  technique.  Assuming  (from  Figure  18)  that  iQ  is  approximately  1000,  this 
would  yield  a  Kolmogorav  length  scale  of  ^lO"^  m.  The  relative  size  of  the  Kolmogorav 
length  scale  to  the  integral  length  scale  is  on  the  order  of  100  times  (the  Kolmogorav 
length  scale  is  100  times  smaller  than  the  integral  length  scale).  This  can  be  checked  by 
using  a  relation  of  the  small  length  scales  (Kolmogorav)  to  large  length  scale  (integral 
length  scale): 


l  =  Rc 

L  ^ 


or 


V 


-1/2 


UA 


X 


For  the  flow  in  this  work,  the  mean  velocity,  U,  is  approximately  20  m/s,  v  is 
approximately  20  xlO"*  m^/s  and  Ax  is  2-3  cm  (depending  on  whether  blowing  is 
initiated).  Using  these  estimated  values,  the  value  of  Re  *^  suggests  that  the  Kolmogorav 
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length  scale  (Ti)  is  approximately  140  times  smaller  than  the  integral  length  scale  (A*), 
which  is  in  good  agreement  with  the  value  determined  using  the  dissipation  wave 
number. 

For  the  data  taken  in  this  work,  a  hotwire  probe  of  length  is  about  1mm.  The 
Kolmogorav  length  scales  we  estimated  to  be  approximately  10^  m,  which  is  10  times 
smaller  than  the  length  of  the  hotwire.  From  this  analysis,  with  the  current  limitations  on 
the  length  of  hotwire  probes,  the  probe  cannot  be  expected  to  measure  length  scales  in  the 
Kolmogorav  scale  range. 


Figure  19.  Con5»arison  of  Dissipation  Spectrum  with  and  without  Grid  Blowing 
Conclusions 

The  results  of  initial  testing  of  an  actively  blown  grid  for  generation  of  large  length-scale, 
high  intensity  freestream  turbulence  have  been  reported.  The  preliminary  testing  of  the 
grid  showed  that  at  the  current  blowing  rates  (blowing  plenum  pressure  of  60  psig),  the 
grid  could  achieve  a  turbulence  intensity,Tu%,  of  17%  and  an  integral  length  scale.  A*,  of 
3.5  cm.  The  desired  intensity  and  length  scale  from  this  grid  are  approximately  20%  and 
4  cm,  respectively.  Continuing  work  will  optimize  the  blowing  rate  and  investigate  the 
uniformity  of  the  flow  in  the  cross-stream  direction  with  blowing.  Also  presented  were 
statistical  and  spectral  analysis  tools  to  fully  characterize  the  turbulence  generated  by  the 
grid,  both  with  and  without  active  blowing.  From  these  analyses,  the  decay  of  turbulence 
intensity  and  length  scale  in  the  streamwise  direction  were  presented  and  compared  to 
published  correlations  for  decay  of  turbulence  and  dispersion  of  length  scale.  The  current 
data  without  blowing  agree  well  with  the  turbulent  decay  correlations,  but  there  is  some 
discrepancy  between  the  measured  length  scale  dispersion  and  the  published  correlations. 
This  discrepancy  will  be  investigated  in  further  work.  The  blown  grid  data  seem  to 
follow  a  different  decay  law,  which  varies  with  blowing  rate. 

Spectral  analysis  of  the  data  demonstrated  that  the  blown  grid  power  spectral  density  and 
energy  content  at  different  wave  numbers  agrees  well  with  the  Von  Karman  spectrum 
equations,  which  are  the  theoretical  energy  spectrum  for  purely  isotropic  turbulence  in 
the  absence  of  outside  forces  and  viscous  dissipation.  The  data  for  the  unblown  grid  was 
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shown  to  deviate  slightly  from  the  Von  Karman  spectrum  in  the  region  prior  to  the 
inertial  subrange.  It  was  shown  that  the  unblown  grid  data  contains  additional  energy  in  a 
small  frequency  range  corresponding  to  a  frequency  where  vortex  shedding  can  be 
expected  from  the  bars.  This  vortex  shedding  was  demonstrated  to  be  relatively  weak 
and  was  shown  to  decay  as  the  flow  progressed  in  the  streamwise  direction.  The  blown 
grid  data  was  also  shown  to  have  a  larger  inertial  subrange  than  the  grid  data  without 
blowing. 

Finally,  an  attempt  was  made  at  determining  the  dissipation  (Kolmogorav)  length  scale. 
Through  analysis  of  the  dissipation  energy  spectra  of  the  data,  the  dissipation  wave 
number  was  estimated.  With  the  data  collected,  it  was  difficult  to  discern  a  difference 
between  the  dissipation  scale  between  the  blown  and  unblown  data.  The  dissq)ative 
length  scales  were  estimated  within  the  limitations  of  the  measurement  technique  to  be  on 
the  order  of  100  times  smaller  than  the  integral  length  scale,  which  is  the  most  energetic 
scale.  ■  . 
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ABSTRACT 

Experiments  were  performed  in  a  transonic  cascade  wind  tunnel 
to  investigate  the  film  effectiveness  and  heat  transfer  coefficient  on  the 
suction  side  of  a  high-turning  turbine  rotor  blade.  The  coolant  scheme 
consisted  of  six  rows  of  staggered,  discrete  cooling  holes  on  and  near 
the  leading  edge  of  the  blade  in  a  showerhead  configuration.  Air  was 
cooled  in  order  to  match  the  density  ratios  found  under  engine 
conditions.  Six  high-frequency  heat  flux  gauges  were  installed 
downstream  of  the  cooling  holes  on  the  suction  side  of  the  blade. 
Experiments  were  performed  with  and  without  film  and  the  coolant  to 
freestream  total  pressure  ratio  was  varied  from  1.02  to  1.19.  In  order 
to  simulate  real  engine  flow  conditions,  the  exit  Mach  number  was  set 
to  1.2  and  the  exit  Reynolds  number  was  set  to  5x10^.  The  freestream 
turbulence  was  approximately  1  %.  The  heat  transfer  coefficient  was 
found  to  increase  with  the  addition  of  film  cooling  an  average  of  14% 
overall  and  to  a  maximum  of  26%  at  the  first  gauge  location.  The 
average  film  cooling  effectiveness  over  the  gauge  locations  was  25%. 
Both  the  heat  transfer  coefficient  and  the  film  cooling  effectiveness 
were  found  to  have  only  a  weak  dependence  upon  the  coolant  to 
freestream  total  pressure  ratio  at  the  gauge  locations  used  in  this  study. 

NOMENCLATURE 

B  blowing  ratio  (pu)c/(pu)f 

d  cooling  hole  diameter  (1  mm  ) 

C  specific  heat  of  air,  1005  J/(kg  K) 

h  heat  transfer  coefficient  w/o  film  cooling 
he  heat  transfer  coefficient  w/  film  cooling 

I  momentum  flux  ratio  (pu^)c/(pu^)f 

M  density  ratio  pc/pf 

p  static  pressure 

Pr  Prandtl  Number,  0.71 

q  heat  flux 


r  recovery  factor 

T.vt,  local  adiabatic  wall  temperature 

Tc  coolant  exit  temperature 

Td  TrTr 

Tf  local  recovery  temperature 

T|  freestream  total  temperature 

Tw  local  wall  or  blade  temperature 

u  local  freestream  velocity 

T|  film  effectiveness 

Subscripts 

c  coolant  or  w/  rilm  cooling 

f  freestream 

INTRODUCTION 

The  efficiency  of  a  gas  turbine  engine  increases  with  turbine  inlet 
temperature.  In  the  ongoing  effort  to  raise  the  turbine  inlet 
temperature  the  gas  stream  temperature  is  made  to  greatly  exceed  the 
operating  temperatures  of  blade  materials,  requiring  elaborate  blade 
cooling  techniques  to  be  developed.  One  such  cooling  scheme  is  film 
cooling,  in  which  cool  air  drawn  from  the  compressor  is  forced 
through  holes  on  the  surface  of  the  turbine  blade  in  order  to  insulate 
the  blade  from  the  hot  engine  environment.  Since  this  air  is  seen  as  a 
loss  of  work  for  the  engine,  it  is  necessary  for  the  thermal  designer  to 
gain  an  understanding  of  the  physics  of  film  cooling  and  the 
parameters  that  affect  it  in  order  to  minimize  the  air  drawn  from  the 
compressor. 

In  an  attempt  to  provide  information  about  the  heat  transfer  to 
film  cooled  turbine  blades,  many  flat  plate  experiments  have  been 
performed.  These  experiments  allowed  researchers  to  amass  a  great 
deal  of  information  about  the  parameters  that  affect  heat  transfer  into  a 
film  cooled  surface.  A  compilation  of  early  experimental  woiic  done 
on  flat  plates  is  given  by  Goldstein  (1971). 


Schwarz  et  al.  (1990)  investigated  the  effects  of  curvature  and 
blowing  ratio  on  the  film  cooling  effectiveness  of  a  convex  surface. 
They  found  that  curved  surfaces  had  better  film  cooling  effectiveness 
than  flat  plates.  Also,  Schwarz  et  al.  found  that  at  low  blowing  rates, 
film  cooling  is  more  effective  on  the  suction  side  of  the  blade  than  on 
the  pressure  side.  This  work  was  an  intermediate  step  between  flat 
plates  and  investigations  of  heat  transfer  into  realistic  blade 
geometries. 

Horton  et  al.  (1985)  investigated  film  cooling  heat  transfer  of  a 
realistic  turbine  blade  profile  mounted  in  a  linear  cascade.  Thin-film 
resistance  sensors  were  used  in  a  short  duration  blow  down  tunnel  to 
capture  the  transient  surface  temperature  history.  Horton  et  al.  were 
able  to  obtain  measurements  of  the  heat  transfer  coefficient  at  several 
locations  along  both  the  pressure  and  suction  sides  of  the  blade. 

Researchers  at  G.  E.  Aircraft  Engines  performed  heat  transfer 
measurements  on  a  heavily  film  cooled  inlet  nozzle  guide  vane  (NGV) 
in  a  linear  cascade  with  realistic  flow  conditions  (Abuaf  et  al.,  1997). 
Thermocouples  embedded  in  a  thin-Walled  turbine  blade  were  used  to 
measure  the  transient  surface  temperature.  Heat  transfer  coefficient 
and  effectiveness  profiles  were  obtained  for  both  sides  of  the  blade. 
Experimental  results  showed  that  the  heat  transfer  coefficient 
increased  as  film  cooling  was  added. 

The  effects  of  freestream  turbulence  on  the  heat  transfer  into  a 
film  cooled  turbine  blade  were  investigated  by  Ekkad  et  al.  (1995). 
Thermochromic  liquid  crystals  (TLC’s)  were  used  to  obtain  the 
transient  surface  temperature  history.  A  general  conclusion  was  shown 
that  high  freestream  turbulence  results  in  an  increased  heat  transfer 
coefficient  and  a  slightly  reduced  film  cooling  effectiveness. 

Researchers  at  Oxford  University  used  thin-film  gauges  to  study 
the  heat  transfer  into  a  heavily  film  cooled  NGV  (Guo  et  al.,  1996  and 
1997).  Their  measurements  were  performed  in  an  annular  cascade 
with  attempts  made  to  simulate  realistic  engine  conditions.  A  mixture 
of  SF^  and  Argon  was  used  as  a  coolant.  The  heat  transfer  coefficient 
and  film  cooling  effectiveness  profiles  were  determined  for  the 
pressure  side  of  the  NGV. 

In  1998,  Drost  et  al.  used  TLC’s  to  record  the  transient  surface 
temperature  history  of  an  NGV  airfoil  exposed  to  a  step  input  in  heat 
transfer  while  in  a  linear  cascade.  They  used  a  foreign  gas  as  a  coolant 
to  obtain  a  realistic  density  ratio.  They  concluded  that  the  film  cooling 
effectiveness  was  higher  near  the  cooling  holes  and  that  the 
mainstream  turbulence  level  had  a  weak  effect  on  suction  side 
effectiveness,  but  that  it  increased  suction  side  heat  transfer 
coefficients. 

The  Virginia  Tech  facility  consists  of  a  transonic  linear  cascade 
of  first  stage  turbine  blades  with  a  shower  head  film  cooling  scheme. 
Preliminary  woric  on  the  same  blade  and  cooling  scheme  with  only  a 
single  gauge  location  was  presented  by  (Popp  et  al.,  1999). 
Measurements  were  made  using  thin-film  heat  flux  gauges  of  both  the 
heat  transfer  coefficient  and  film  cooling  effectiveness  at  six  locations 
along  the  blade  suction  side  surface.  In  the  present  paper, 
measurements  along  the  length  of  the  chord  were  made,  and  the 
influence  of  blowing  parameters  on  the  heat  transfer  coefficient  and 
film  cooling  effectiveness  is  presented.  Part  II  of  this  research 
presents  the  unsteady  effects  of  shocks. 

EXPERIMENTAL  APPARATUS 

Transonic  Wind  Tunnel 

The  experiments  described  in  this  paper  were  all  performed  in  the 
Virginia  Tech  transonic  wind  tunnel,  an  intermittent  blow  down 
facility  with  an  open  discharge.  A  reciprocating  compressor  is  used  to 
load  a  high-pressure  tank  that  serves  as  the  supply  for  the  mainstream 


flow.  When  released,  air  flows  at  a  rate  of  approximately  10  kg/s 
through  the  system  of  pipes  shown  in  Fig.  1.  A  fast-acting  control 
valve  is  used  to  maintain  the  tunners  total  pressure  at  a  constant  value 
that  was  experimentally  determined  to  give  the  desired  exit  Mach 
number.  The  pressure  is  constantly  monitored  and  controlled  by  a 
pressure  transducer  upstream  of  the  test  section.  Tbe  Virginia  Tech 
wind  tunnel  also  has  the  capacity  for  heated  flows  by  utilizing  the 
lower  loop  shown  in  Fig.  1.  In  the  bottom  of  the  loop  there  are  two  36 
kW  resistance  heaters  which  heat  the  air.  This  heated  air  is  then 
circulated  over  a  bank  of  copper  tubes  in  the  top  of  the  loop  by  a  small 
fan.  The  heated  copper  tubes  act  as  a  passive  heat  exchan:^er  and  heat 
the  tunnel* s  mainstream  flow  during  a  run.  During  the  course  of  a 
tunnel  run  the  mainstream  flow  temperature  drops,  as  will  be 
described  later,  this  allows  a  range  of  density  ratios  over  which  the 


Figure  2:  Schematic  of  Cascade  Test  Section 


Cascade  Test  Section 

The  test  section  used  to  mount  the  rotor  blades  was  constructed  of 
aluminum  with  clear-acrylic  side  walls.  The  side  walls  give  optical 
access  to  the  flow  in  the  passages  and  allow  shadowgraph  and 
Schlieren  photography  to  be  performed.  A  cut-a-way  of  the  test 
section  can  be  seen  in  Fig.  2.  TTie  test  section  contains  four  full  blades 
and  two  half  blades,  a  total  of  five  flow  passages.  Flow  enters  the  test 
section  at  Mach  0.3,  is  turned  128®  while  passing  over  the  blades,  and 
exits  the  test  section  at  Mach  1.2.  The  shaded  blade  was  instrumented 
and  used  to  acquire  all  of  the  data  presented  in  this  paper.  Schlieren 
imaging  was  used  to  demonstrate  flow  periodicity  in  the  two  passages 
around  the  instrumented  blade. 
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Turbine  Blade  Design 

The  blade  profile  used  in  this  study  was  a  generic  design  given  by 
GE  Aircraft  Engines  and  was  made  from  aluminum.  Measured  and 
predicted  Mach  number  profiles  matched  well  for  this  blade  design 
and  were  presented  by  Hale  et  al.  (1997).  The  Reynolds  number  based 
on  aerodynamic  chord  (13.6  cm)  and  unheated  exit  conditions  is  about 
6x10^.  The  maximum  value  of  the  acceleration  factor  in  the  region  of 
the  sensors  was  approximately  l.lxlO"^.  Six  rows  of  cooling  holes  are 
located  near  the  leading  edge  in  a  showerhead  arrangement,  as  can  be 
seen  in  Fig.  3.  All  of  the  holes  have  a  diameter  of  1  mm  and  the 
spacing  Ltitween  holes  is  approximately  9. 1  hole  diameters.  The  four 
rows  of  cooling  holes  nearest  the  leading  edge  are  all  aligned  normal 
to  the  streamwise  tangent  and  are  inclined  at  an  angle  of  30®  to  the 
surface  where  they  emerge.  The  two  rows  of  coolant  holes  farthest 
from  the  leading  edge  are  both  aligned  parallel  to  the  streamwise 
tangent.  The  gill  holes  on  the  suction  side  are  inclined  at  an  angle  of 
30®  to  the  surface  where  they  emerge,  while  the  gill  holes  on  the 
pressure  side  are  inclined  at  an  angle  of  45®  to  the  surface  where  they 
emerge. 
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Figure  3:  Blade  Cooling  Scheme 


Senso 

Six  icis  of  gauges  were  press-fit  into  the  blade  downstream  of  the 
last  row  of  cooling  holes  on  the  suction  side  as  can  be  seen  in  Fig.  4. 
Figure  4  is  a  plan  view  of  the  suction  side  of  the  blade  showing  the 
locations  of  the  cooling  hole  exits  and  gauge  locations.  Each  set  of 
gauges  is  composed  of  a  high-speed  pressure  transducer  (Kulite  XCQ- 
062-50a),  a  high-speed  heat  flux  and  temperature  sensor  (HFM-7),  and 
a  surface-mounted  thermocouple.  The  Kulite  transducers  are  used  to 
acquire  high-speed  static  pressure  measurements  on  the  blade  surface. 
The  HFM-7 ’s  acquire  separate,  direct  and  simultaneous  measurements 
of  both  the  local  heat  flux  and  wall  temperature.  The  HFM-7  has  a 
gauge  diameter  of  5mm  and  an  active  diameter  of  4mm.  A  detailed 
discussion  of  these  gauges  can  be  found  in  Diller  (1993),  The  surface- 
mounted  thermocouple  is  used  as  a  calibration  tool  for  the  HFM-7’s 
temperature  sensor.  The  first  set  of  gauges  was  mounted 
approximately  9.1  hole  diameters  downstream  of  the  last  row  of 
cooling  holes  and  the  subsequent  sets  of  gauges  were  staggered  along 
the  span  and  in  the  streamwise  direction.  All  of  the  gauges  were 
mounted  within  a  three- inch  wide  section  of  the  blade  span  which 
surface  oil  flow  visualizations  had  shown  to  be  free  of  endwall  effects. 
The  high  thermal  conductivity  gauge  substrates  (aluminum  nitride) 
were  press  fit  into  the  aluminum  blade  to  minimize  the  fluid  and 
therm^  disruptions  of  the  surface.  The  effect  of  a  possible  sensor 
mismatch  with  the  surface  contour  was  investigated  by  Peabody  and 


Diller  (1998)  and  deemed  to  be  negligible  for  this  case.  Coolant 
temperature  measurements  were  made  with  very  small  exposed 
junction  thermocouples  protruding  into  the  exit  of  the  last  set  of 
cooling  holes.  The  values  were  mass-averaged  and  compared  with  the 
coolant  temperature  in  the  insulated  plenum.  The  temperatures 
throughout  the  coolant  supply  were  tested  for  uniformity. 
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Figure  4:  Plan  View  of  the  Suction  Side  Showing  Blade 
Instrumentation  and  (I!oolant  Exits 

Coolant  Supply  System 

A  schematic  of  the  coolant  system  is  shown  in  Fig.  5.  Air  is 
supplied  to  the  insulated  plenum  of  the  instrumented  blade  from  a 
large  storage  tank.  The  coolant  supply  storage  tank  is  fed  by  a  5  hp 
compressor  with  a  filter  and  dryer  system  attached.  When  air  is 
released  from  the  storage  tank  it  passes  through  a  coolant  control 
system.  This  coolant  control  maintains  the  plenum-to-freestream  total 
pressure  ratio  over  the  course  of  the  run  via  an  integral  feedback 
control.  Controlling  the  pressure  ratio  allows  the  momentum  ratio  to 
be  controlled,  but  the  density  and  blowing  ratios  vary  over  a  run.  The 
air  then  flows  through  an  orifice  plate,  where  the  mass  flow  rate  is 
monitored,  and  into  a  chiller  which  uses  liquid  nitrogen  as  a  coolant. 
The  cold  air  is  then  fed  into  the  insulated  plenum  of  the  blade. 


Figure  5:  Schematic  of  Coolant  System 


DATA  ANALYSIS  AND  RESULTS 

Three  different  types  of  experiments  are  presented  here.  The  first 
type,  uncooled  runs,  used  a  plenum  that  was  plugged  with  a  plastic 
insert  to  block  coolant  flow  during  the  run  and  to  insure  that  no 
mainstream  flow  bled  through  different  rows  of  coolant  holes.  The 
second  type,  cooled  runs,  used  the  coolant  flow  at  the  design  pressure 
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ratio  of  1 .04.  Lastly,  a  third  set  of  experiments  were  performed  where 
the  pressure  ratio  was  varied  over  the  range  1 .02  to  1 .20  in  order  to 
investigate  effect  of  pressure  ratio  on  both  the  heat  transfer  coefficient 
and  film  cooling  effectiveness. 

Uncooied  Data  Analysis  and  Results 

The  dominant  mode  of  heat  transfer  in  gas  turbines  is  convection, 
and  the  convection  heat  transfer  coefficient  is  defined  by  the  equation, 

q  =  h(T.,-Tj  (1) 

where  the  difference  between  the  adiab^it'C  wall  temperature,  T»w»  and 
the  wall  temperature,  Tw,  is  the  driving  force  of  the  heat  transfer.  For 
high-speed  flows  with  no  film  inj^^uon,  is  equivalent  to  the 
recovery  temperature,  T^,  and  equation  (1)  can  be  rewritten  as, 

q  =  h(T,-T,)  (2) 

The  recovery  temperature  can  be  related  to  the  freestream  total 
temperature  by  the  recovery  factor  and  the  freestream  velocity, 

2 

T,-T,  =(l-r)~  =  T,  (3) 

2C 

where  r  is  the  recovery  factor,  u  is  the  freestream  velocity,  and  C  is  the 
specific  heat  of  air.  The  difference  between  the  freestream  total 
temperature  and  the  recovery  temperature  is  a  constant  during  the  run 
and  will  hereafter  be  referred  to  as  Td.  The  recovery  factors  for 
laminar  and  turbulent  boundaiy  layers  are  usually  considered  to  be 
Pr*^  and  Pr*^,  respectively  [Kays  and  Crawford,  1993].  Substituting 
equation  (3)  into  (2)  yields, 

q  =  h(T,-T„)-hT,  (4) 

In  equation  (4),  the  local  heat  flux,  q,  is  the  dependent  variable 
and  the  temperature  difference,  -  T^,  is  the  independent  variable. 
The  heat  transfer  coefficient  is  the  slope  in  this  equation  and  Td  is 
equal  to  the  temperature  difference,  T^  -  T^,  when  q  =  0  (i.e.  the  x-axis 
intercept).  It  was  shown  by  Popp  et  al.  (1999)  that  the  heat  transfer 
coefficient  and  Td  did  not  significantly  depend  upon  the  temperatures 
involved  in  equation  (4).  Therefore,  for  the  uncooled  case,  the  heat 
transfer  coefficient  and  the  temperature  difference  Td  can  be  obtained 
by  measuring  the  freestream’ s  total  tei,iperature,  the  blade’s  wall 
temperature,  and  the  local  heat  flux  ii.to  the  blade  over  the  course  of 
an  experiment. 


Figure  6:  Uncooled  Data  Analysis  (Gauge  Location  1) 


When  the  experimental  data  was  fitted  to  equation  (4),  the  result 
was  a  linear  trend  as  can  be  seen  in  Fig.  6.  Figure  6  is  a  typical  result 
of  the  analysis  performed  on  data  gathered  at  the  first  gauge  location. 
A  regression  analysis  was  performed  to  fit  a  straight  line  through  this 
data.  The  heat  transfer  coefficient  was  obtained  from  the  slope  of  this 
line,  and  the  temperature  difference  Td  was  determined  from  the  x-axis 
intercept. 

A  typical  time  history  of  all  relevant  temperatures,  heat  flux  and 
the  heat  transfer  coefficient  from  gauge  location  1  can  be  seen  in  Fig. 
7.  The  recovery  temperature  in  this  figure  is  determined  as  a  drop 
from  the  measured  freestream  temperature  using  the  calculated  local 
value  of  Td  (shown  in  Fig  6).  The  freestream’ s  total  temperature  is 
seen  to  vary  significantly  over  the  course  of  the  experiment. 
Typically,  it  varies  from  100®C  to  30®C  as  the  passive  heat  exchanger 
cools  down.  The  blade  wall  temperature  rises  ^proximately  8®C  over 
the  course  of  the  run.  The  time  history  of  the  heat  transfer  coefficient 
was  calculated  from  the  corresponding  time  histories  of  Tr  T^,  and  q, 
using  equation  (2).  It  can  be  seen  that  the  heat  flux  into  the  blade 
approaches  zero  as  the  recovery  and  wall  temperatures  converge, 
which  is  consistent  with  the  definition  of  heat  flux  in  equation  (2). 


Figure  7:  Time  History  of  T| ,  T^ ,  T^ ,  q,  and  h 
(Uncooled  Run,  Dashed  Lines  Indicate  Useful  Range  of  Data) 

Six  experimental  runs  were  performed  without  film  cooling. 
Figure  8  shows  the  heat  transfer  coefficients  determined  at  each  gauge 
location  for  all  six  runs.  The  average  level  of  the  heat  transfer 
coefficient  over  the  region  of  interest  was  700  W/m^-®C.  The  trend  of 
the  heat  transfer  coefficient  profile  was  suspected  to  be  caused  by 
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transition  of  the  boundary  layer.  However,  the  same  trend  was  still 
observed  during  experiments  with  the  suction  side  boundary  layer 
tripped  at  the  leading  edge. 


Figure  8:  Heat  Transfer  Coefficient  Profile  (Uncooled  Run) 

Film  Cooled  Data  Analysis  and  Results 

For  the  film  cooled  experiments,  it  was  important  to  be  able  to 
visualize  the  condition  of  the  coolant  film.  The  optical  access  of  the 
cascade  test  section  allowed  both  shadowgraph  and  Schlieren 
photographs  to  be  made  in  order  to  visualize  the  film  cooling  layer.  A 
high-speed  digital  camera  was  used  to  capture  the  shadowgraph  shown 
in  Fig.  9.  The  coolant  film  is  clearly  attached  and  turbulence  is 
evident  in  the  film  layer  on  the  suction  side.  The  pressure  side  of  the 
blade  is  obscured  by  equipment  outside  of  the  test  section.  Although 
not  evident  in  this  photograph,  the  stagnation  point  has  shifted  slightly 
towards  the  suction  side,  causing  the  coolant  jets  exiting  the  nose  holes 
to  be  diverted  entirely  to  the  pressure  side.  This  did  not  noticeably 
alter  the  pressure  distribution  around  the  blade,  but  none  of  the  coolant 
from  the  stagnation  holes  flowed  to  the  suction  side  of  the  blade  for 
cooling.  Therefore,  only  the  three  rows  of  coolant  holes  on  the  suction 
side  provide  protection  for  the  suction  side  of  the  blade. 


Figure  9:  Shadowgraph  of  Cooling  Film 


With  the  addition  of  film  cooling  to  the  surface  of  a  turbine  blade 
in  a  high-speed  flow,  the  governing  equation  of  convection  heat 
transfer  is  once  again  equation  (1).  However,  in  such  cases  the 
adiabatic  wall  temperature  is  usually  non-dimensionalized  using  the 
adiabatic  film  effectiveness  defined  as 


where  r\  is  the  adiabatic  film  effectiveness.  Equation  (5)  can  be  solved 
for  T,w  and  then  substituted  into  equation  (1)  which  will  yield, 


The  dependent  variable  in  equation  (6)  is  the  heat  flux  divided  by  a 
temperature  difference  and  the  independent  variable  is  a  ratio  of  two 
temperature  differences.  When  the  experimental  results  are  plotted 
using  these  variables,  the  slope  of  equation  (6)  is  the  heat  transfer 
coefficient  and  the  x-axis  intercept  gives  the  adiabatic  film 
effectiveness.  Using  the  values  of  T^  determined  for  each  gauge 
position  in  the  uncooled  case,  the  recovery  temperatures  at  each 
location  were  calculated  from  the  freestream  total  temperature. 

The  experimental  data  from  a  typical  film  cooled  run  were  fitted 
to  equation  (6)  and  the  result  is  a  linear  trend,  which  is  displayed  in 
Fig.  10.  A  regression  analysis  was  performed  to  fit  a  line  through  the 
data.  The  film  cooled  heat  transfer  coefficient  was  determined  from 
the  slope  of  the  line  and  the  adiabatic  film  effectiveness  was  obtained 
from  the  x-axis  intercept.  The  resulting  values  do  not  appear  to  be  a 
function  of  temperature.  For  the  same  gauge  location  as  shown  earlier 
in  the  uncooled  case,  the  heat  transfer  coefficient  has  increased  to  950 
W/m^  ®C,  while  the  effectiveness  is  31%. 


Figure  10:  Film  cooled  Data  Analysis 

A  sample  time  history  of  a  film  cooled  run  can  be  seen  in  Fig.  11. 
The  test  window  indicates  the  limits  of  the  data  set  that  were  used  in 
the  film  cooled  analysis.  The  traces  in  Fig.  1  la  are  similar  to  those  in 
Fig.  7a  except  for  the  addition  of  a  mass-averaged  coolant 
temperature.  The  adiabatic  wall  temperature  is  calculated  using  the 
film  effectiveness  determined  from  Fig.  10  and  equation  (5).  The 
freestream  temperature  again  falls  from  around  100*^  to  30®C,  while 
the  blade  temperature  rises  from  near  0®C  to  about  10®C.  The  liquid 
nitrogen  heal  exchanger  allows  the  coolant  temperature  to  be 
maintained  at  temperatures  ranging  between  -120®C  and  -150°C.  In 
comparing  Figs.  11a  and  7a,  it  should  be  noted  that  the  driving 
temperature  of  the  heat  transfer  is  lower  for  the  film  cooled  case  than 
for  the  uncooled  case.  Figure  11b  shows  that  even  though  the  heat 
transfer  coefficient  at  location  1  increases,  the  adiabatic  wall 
temperature  is  lowered  enough  that  the  peak  level  of  the  heat  flux 
decreases  from  4.5  W/cm^  to  3.5  W/cm^. 
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Figure  1 1 :  Time  Histories  for  a  Film  cooled  Run  (a)  All  Relevant 
Temperatures 

(b)  Heat  Flux  and  Heat  Transfer  Coefficient 
(c)  Density  Ratio,  Blowing  Ratio,  and  Momentum  Ratio 
(Dashed  Lines  Indicate  Useful  Range  of  Data) 


the  film  cooled  heat  transfer  coefficient  seen  in  Fig.  1  lb.  A  singularity 
occurs  in  the  calculation  of  the  film  cooled  heat  transfer  coefficient 
when  T,w  and  Tw  converge.  Aside  from  the  singularity,  the  film  cooled 
heat  transfer  coefficient  is  fairly  constant  over  the  course  of  the 
experiment. 

Figure  11(c)  shows  the  coolant-to-fi-eestream  density  ratio  and 
pressure  ratio.  In  addition,  Fig.  1 1(c)  shows  the  blowing  ratio  and  the 
momentum  ratio  for  the  three  coolant  holes  that  affect  the  suction  side 
of  the  blade.  All  of  these  ratios  were  determined  assuming  isentropic 
flow  conditions  in  the  freestream. 

Six  film  cooled  exr'i  rimental  runs  were  performed  and  analyzed 
in  the  manner  just  described.  Figure  12  shows  the  results  of  heat 
transfer  coefficient  and  film  effectiveness  for  all  six  runs.  The  average 
level  of  the  heat  transfer  coefficient  for  the  film  cooled  case  was  found 
to  be  850  W/m^  ®C,  a  21%  increase  over  the  uncooled  experiments. 
The  largest  increase  was  at  gauge  location  1  followed  by  2  and  6.  The 
smallest  increase  due  to  film  cooling  was  at  gauge  location  3.  Over  the 
region  studied,  the  film  had  an  average  adiabatic  effectiveness  of  25%. 
Possible  explanations  (e.g.,  transition,  gauge  calibration,  span  wise 
variations)  for  the  minimum  h  and  maximum  T|  at  gauge  location  3 
were  checked,  but  none  appeared  to  be  the  cause.  These  experiments 
provide  a  basis  for  the  unsteady  work  that  will  be  investigated  in  Part  2 
(Popp,  2000). 


The  adiabatic  wall  temperature  was  used  in  equation  (1)  along 
with  the  heat  flux  and  wall  temperature  to  generate  the  time  history  of 


Figure  12:  Film  Cooled  Experiments  at  a  Pressure  Ratio  of  1.04 
(a)  Heat  Transfer  Coefficient,  (b)  Adiabatic  Film  Effectiveness 
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UNCERTAINTY  IN  h  AND  Ti 

The  theoretical  uncertainty  for  the  h  and  Tj  measurements  are 
divided  into  two  categories;  bias  and  precision  uncertainties.  The  bias 
uncertainty  is  composed  entirely  of  the  calibration  of  the  heat  flux 
microsensors.  The  uncertainties  and  calibration  techniques  for  these 
gauges  are  detailed  in  a  paper  by  Smith  et.  al.  (1999).  Separate 
calibrations  were  done  using  convection  and  radiation  and  a  transient 
technique  was  used  to  check  the  calibration  of  each  gauge  at  the  start 
of  every  test.  Table  1  details  these  bias  uncertainties  (based  on  a  90% 
confidence  interval  of  a  student-T  distribution). 


Gauge 

Bias  Uncertainty 

1 

17.0% 

2 

3.7% 

3 

5.2% 

4 

15.2% 

5 

8.8  % 

6 

10.7% 

Table  1 :  Heat  Transfer  Coefficient  Bias  Uncertainties 


These  gauge  bias  uncertainties  do  not  effect  the  bias  uncertainty  in  T| 
since  Ti  is  determined  solely  from  temperature  measurements  (i.e.  the 
x-axis  intercept  in  Fig.  10). 

Precision  uncertainties  are  comprised  of  the  theoretical  variations 
of  temperature  measurements  and  a  scheme  outlined  by  Moffat  (1988) 
for  calculations  involving  complex  computer  algorithms.  Based  on  an 
individual  temperature  uncertainty  of  1.1  °C,  this  type  of  analysis 
yields  the  precision  uncertainties  detailed  in  Table  2. 


Gauge 

h  Precision 
Uncertainty 

T|  Precision 
Uncertainty 

1 

5.9% 

18.5  % 

2 

6.1  % 

21.0% 

3 

5.6% 

13.9% 

4 

5.9% 

21.3% 

5 

6.0% 

17.0% 

6 

5.7% 

11.5% 

Table  2:  Precision  Uncertainties 


Another  indication  of  measurement  uncertainty  is  an 
experimental  value  of  measurement  repeatability.  This  repeatability 
value  was  determined  from  a  number  of  different  experiments 
performed  at  the  same  test  conditions.  The  results  shown  in  Table  3 
are  based  on  8  runs  at  a  pressure  ratio  of  1 .04  and  a  90%  confidence 
bounds  using  a  student-T  distribution. 


Gauge 

h  Repeatability 

Ti  Repeatability 

1 

11.1  % 

16.4% 

2 

4.0% 

8.2% 

3 

5.5% 

7.1  % 

4 

5.1  % 

15.8% 

5 

6.7% 

21.1  % 

6 

7.3% 

11.4% 

Table  3:  Measurement  Repeatability 


THE  INFLUENCE  OF  PRESSURE  RATIO  ON  h  AND  T] 

Experiments  were  performed  using  the  above  data  reduction 
technique  to  determine  the  influence  of  the  plenum  to  freestream  total 
pressure  ratio  on  both  h  and  T|.  Using  the  coolant  control  feedback 
system,  the  pressure  ratio  was  varied  from  run  to  run  over  the  ranges 


of  1.02  to  1.19.  The  results  from  a  number  of  runs  are  shown  in  Fig. 
13  and  Fig,  14,  The  error  bars  on  these  plots  indicate  the  measurement 
repeatablilty  shown  in  Table  3. 


Pressure  Ratio  [  -  ]  Pressure  Ratio  [  -  ] 

Figure  13:  The  Influence  of  Pressure  Ratio  on  the  Heat  Transfer 
Coefficient 


Pressure  Ratio  [  -  ]  Pressure  Ratio  [  -  ] 


1.05  1.1  1.15 

Pressure  Ratio  [  ] 

Figure  14:  The  Influence  of  Pressure  Ratio  on  Film  Cooling 
Effectiveness 
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From  Fig.  13  it  can  be  seen  that,  in  general,  the  heat  transfer 
coefficent  is  fairly  independent  of  the  pressure  ratio.  Also  Fig.  14 
shows  an  independence  of  pressure  ratio  and  effectiveness  for  gauge 
locations  2,  4,  5,  and  6  and  a  slight  decrease  of  effectiveness  at 
locations  1  and  3.  The  slight  decrease  at  these  locations  is  suspected 
to  be  attributed  to  spanwise  shifting  of  the  coolant  at  higher  pressure 
ratios  which  yielded  lower  film  coverage  at  these  locations. 

CONCLUSIONS 

Heat  transfer  experiments  were  performed  on  the  suction  surface 
of  a  fl.n  cooled  transonic  turbine  rotor  blade.  A  method  that  allowed 
for  the  determination  of  both  the  heat  transfer  coefficient  and  either  the 
recovery  temperature  (in  uncooled  runs)  or  the  film  cooling 
effectiveness  (film  cooled  runs)  from  one  transient  run  was  used.  For 
the  nominal  pressure  ratio  of  1.04,  an  average  increase  in  the  heat 
transfer  coefficient  of  2 1  %  was  seen  with  the  addition  of  film  cooling. 
Varying  the  pressure  ratio  in  the  range  from  1.02  to  1.19  had  only  a 
weak  effect  on  the  heat  transfer  coefficient  and  slightly  decreased  the 
film  effectiveness  at  some  gauge  locations. 
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ABSTRACT 

This  paper  reports  on  an  investigation  of  the  heat  transfer  on  the 
suction  side  of  a  transonic  film  cooled  turbine  rotor  blade  in  a  linear 
cascade.  A  shock  wave  is  generated  in  a  shock  tube  and  passes  into  the 
cascade  upstream  of  the  blade  row  to  simulate  the  interaction  of  a 
shock  emerging  ftom  the  trailing  edge  of  an  upstream  nozzle  guide 
vane  on  the  downstream  rotating  blade  row.  The  unsteady  heat  transfer 
caused  by  the  shock  wave  passing  over  the  suction  side  of  the  rotor 
blade  is  analyzed  and  compared  to  predictive  models  published  in  the 
literature.  The  investigation  includes  experiments  with  and  without 
film  cooling  and  a  variation  of  the  most  relevant  physical  parameters. 
From  the  analysis  and  the  comparison  r,f  the  unsteady  heat  flux  it  is 
concluded  that  the  heat  transfer  coericieht  is  not  affected  strongly  by 
the  passing  shock.  Most  of  the  unstea-i^  heat  transfer  is  induced  by  the 
variation  of  temperature  caused  by  the  shock  wave. 

NOMENCLATURE 

Cp  specific  heat  of  air 

h  heat  transfer  coefficient  w/o  film  cooling 
he  heat  transfer  coefficient  w/  film  cooling 

k  thermal  conductivity  of  air 

m  counter  in  Eq.  9 

n  counter  in  Eq.  9 

p  static  pressure 

q  heat  flux  per  unit  area 

qc  unsteady  component  of  heat  flux  (in  Rigby’s  Model) 

qh  unsteady  component  of  heat  flux  (in  Rigby’s  Model) 

qM  unsteady  component  of  heat  flux  (in  Moss’  Model) 

T,^  local  adiabatic  wall  temperature 

Tc  coolant  exit  temperature 

Th  unsteady  component  of  temperature  (in  Rigby’s  Model) 

Tm  unsteady  component  of  temperature  (in  Moss*  Model) 

Tr  local  recovery  temperature 

Tc  freestream  total  temperature 

Tw  local  wall  or  blade  temperature 


At  time  between  data  points 

T|  film  effectiveness  defined  in  Eq.(7) 

Superscripts 

"  unsteady  variation 

INTRODUCTION 

The  quest  for  higher  thrust  to  weight  ratios  in  the  development  of 
aero-engines  has  led  to  the  design  of  nozzle  guide  vanes  (NGV)  with 
supersonic  exit  velocities  and  elaborate  film  cooling  schemes. 
Consequently,  the  rotor  blades  are  not  only  subject  to  wake  but  also 
shock  impingement  as  they  pass  behind  the  NGV’s  at  very  high 
rotational  speed.  The  effect  of  this  unsteady  process  on  the  heat 
transfer  to  the  rotor  blade  has  been  the  topic  of  a  variety  of  research 
programs. 

Researchers  at  Texas  A&M  (Ou  et  al.  (1994),  Mehendale  et  al. 
(1994),  Du  et  al.  (1997))  published  time  averaged  results  from  a  low 
speed  cascade  with  and  without  film  cooling.  Wakes  were  produced 
using  a  spoked  wheel  rotating  upstream  of  the  linear  cascade.  It  was 
concluded  that  the  presence  of  the  wakes  increases  the  heat  transfer 
coefficient  only  on  the  downstream  suction  side  by  affecting  the 
transitional  behavior  of  the  boundary  layer.  With  the  introduction  of 
film  cooling  on  the  suction  side,  the  boundary  layer  is  tripped  early, 
and  the  effect  of  the  wakes  decreases  accordingly. 

Heidmann  et  al.  (1997)  presented  experimental  woric  on  a 
stationary  annular  cascade  consisting  of  rotor  blades  with  an  inlet 
Mach  Number  of  0.27.  Wakes  were  generated  using  a  rotating  bar 
mechanism  rotating  in  front  of  the  cascade.  Time-averaged  film 
cooling  effectiveness  was  measured  for  a  showerhead  film  cooling 
scheme  with  two  rows  of  holes  at  the  suction  side  and  the  pressure  side 
and  one  row  in  the  geometric  stagnation  point.  The  wake-induced 
mixing  decreased  the  film  cooling  effectiveness  all  along  the  blade. 
The  effect  was  nK)re  pronounced  on  the  suction  side  since  the  swirl 
generated  by  the  rotating  rods  shifted  the  stagnation  point.  Therefore, 
more  coolant  was  available  on  the  pressure  side,  almost  outweighing 
the  decrease  in  film  effectiveness  due  to  the  increased  mixing. 


TTie  last  in  a  series  of  publications  on  heat  flux  and  pressure 
measurements  in  a  full  stage  rotating  turbine  done  at  the  Calspan-UB 
Research  Center  featured  phase-locked  measurements  of  heat  flux  on 
the  rotor  blades  (Dunn  (1989)).  No  on-rotor  pressure  measurements 
were  provided  to  relate  the  heat  flux  to  the  flowfield.  On  the  shroud 
both  pressure  and  heat  flux  data  were  recorded.  The  fluctuations  of 
heat  flux  and  pressure  were  shown  to  be  in  phase  and  of  similar  shape. 

In  Guenette  et  al.  (1989)  a  first  set  of  time-resolved  on-rotor  heat 
flux  and  pressure  data  from  a  full  stage  blowdown  turbine  was 
reported.  The  mid  section  of  the  rotor  blade  profile  was  identical  to  the 
blade  tested  in  a  linear  cascade  at  the  University  of  Oxford.  The  da^a 
from  both  setups  were  compared  and  showed  reasonable  agreement.  In 
1996  Abhari  et  al.  (1996)  compared  the  heat  flux  results  from  a  ste??dy- 
state  code  to  time-averaged  measurements  on  a  film  cooled  blade.  The 
presence  of  the  wakes  and  shocks  seemed  to  decrease  the  heat  transfer 
on  the  suction  side  and  increase  it  on  the  pressure  side. 

The  most  comprehensive  body  of  research  was  published  by  the 
University  of  Oxford.  In  1985  Doorly  et  al.(1985)  were  the  first  to 
report  on  time-resolved  heat  transfer  measurements  in  a  linear  cascade. 
A  combination  of  wakes  and  weak  shocks  was  produced  using  a 
rotating  bar  mechanism.  Turbulent  spots  and  bubbles  created  by 
shocks  and  wakes  were  shown  to  be  responsible  for  an  increase  of  heat 
flux.  Johnson  et  al.  (1988)  used  slightly  modified  blades  and  increased 
the  rotational  speed  of  the  bars  to  create  stronger  shock  waves.  The 
modified  blades  were  the  same  ones  used  by  Guenette  et  al.  (1989). 
The  unsteady  heal  flux  caused  directly  by  the  impact  of  the  shock 
wave  was  observed  to  be  the  major  contribution  to  the  overall 
unsteady  heat  flux.  A  first  version  of  a  predictive  model  linking 
isentropic  shock  heating  and  one  dimensional  conduction  in  the  fluid 
to  the  unsteady  heat  flux  showed  fairly  good  agreement  with  measured 
heat  flux  data.  Rigby  et  al.  (1989)  modified  this  model  and  gave  a 
thorough  derivation.  Johnson  et  al.  (1990)  found  that  while  the 
prediction  yielded  good  agreement  on  the  suction  side,  the  pressure 
side  unsteady  heat  flux  was  undeipredicted.  The  reason  for  this  was 
found  in  a  “Vortical  Bubble”  created  by  the  shock  wave  and  traveling 
down  the  pressure  side. 

Rigby  et  al.  (1990)  first  included  film  cooling  in  the  Oxford 
cascade  setup  for  the  investigation  of  unsteady  heat  flux.  The  shock 
wave  signature  was  found  to  change  with  the  introduction  of  fil  i 
cooling.  Data  recorded  in  an  entirely  new  full  stage  rotating  turbine 
was  first  published  in  1995  by  Moss  et  al.  (1995).  Time-resolved  on- 
rotor  heat  flux  and  pressure  measurements  were  presented.  It  was 
shovm  that  the  unste^y  heat  flux  could  be  predicted  by  multiplying  a 
constant  heat  transfer  coefficient  and  a  temperature  variation 
calculated  from  the  unsteady  pressure  based  on  an  isentropic 
assumption.  Moss  et  al.  (1997)  compared  time-averaged  heat  flux 
measurements  on  the  rotor  with  the  presence  of  upstream  wakes  and  a 
flowfield  undisturbed  by  the  nozzle  guide  vanes.  The  values  of  time- 
averaged  heat  transfer  with  and  without  the  presence  of  upstream 
disturt)ances  were  nearly  identical,  thus  proving  the  assumption  of  a 
constant  heat  transfer  coefficient  valid. 

Researchers  at  Virginia  Tech  have  studied  in  a  stationary  cascade 
the  individual  components  of  the  disturbance  from  upstream  rotating 
blades.  The  effects  of  velocity  defect,  freestream  turbulence,  and 
transient  shock  impingement  were  measured  in  separate  experiments 
to  understand  the  effect  of  each  on  turbine  blade  heat  transfer.  Hale  et 
al.  (1997)  modeled  the  effect  of  wake  passing  in  a  quasi-steady  way 
using  a  stationary  strut  that  added  turbulence  and  velocity  defect. 
Increases  in  heat  transfer  coefficient  were  measured  for  a  number  of 
locations  on  the  blade,  particularly  on  the  pressure  side.  Nix  et  al. 
(1997)  analyzed  in  detail  the  progression  of  a  shock  wave  through  the 
same  cascade  and  its  effect  on  the  unsteady  heat  transfer.  When 


averaged  over  a  200  ps  blade  passing  event,  a  maximum  increase  of 
heat  flux  of  60%  was  measured  due  to  shock  passing.  Popp  et  al. 
(1999)  presented  time-resolved  heat  flux  and  pressure  data  from  one 
measurement  location  on  the  suction  side  of  a  film  cooled  transonic 
rotor  cascade.  It  was  shown  that  the  heat  transfer  coefficient  as  well  as 
the  film  cooling  effectiveness  were  not  affected  significantly  by  the 
passing  shock  wave. 

The  intent  in  the  present  study  is  to  extend  the  investigation  of 
shock  effects  to  a  larger  number  of  measurement  locations  and  to 
verify  the  conclusions  for  a  variety  of  physical  parameters  for  a  film 
cooled  blade.  The  unsteady  heat  flux  will  be  analyzed  to  find  the 
relative  magnitudes  of  unsteady  heat  transfer  coefficient,  unsteady 
temperature  and  film  effectiveness  and  their  contributions  to  the 
overall  unsteady  heat  flux.  Also,  the  experimental  data  will  be 
compared  to  predictive  models  developed  by  Rigby  et  al.  (1989)  and 
Moss  et  al.  (1995). 

EXPERIMENTAL  APPARATUS 


Wind  Tunnel  Facliitv,  Cascade  and  Shock  Apparatus 

The  experiments  necessary  for  this  investigation  were  performed 
in  the  transonic  blowdown  wind  tunnel  at  Virginia  Tech.  A  passive 
heating  device  is  available  to  achieve  high  (120°C)  inlet  temperatures 
to  the  cascade.  It  consists  of  many  copper  tubes  that  are  preheated 
prior  to  running  the  tunnel.  Fig.  1  of  Part  I  (Smith  et  al.  (2000))  shows 
the  wind  tunnel  with  the  heating  loop.  With  the  present  cascade,  the 
facility  allows  run  times  of  up  to  35  seconds  with  the  inlet  pressure 
controlled.  The  test-section  and  cascade  built  for  this  investigation  are 
shown  in  Fig.  2  of  Part  I  (Smith  et  al.  (2000)).  The  cascade  consists  of 
four  full  and  two  half  blades  forming  five  passages.  The  blade  design 
is  a  generic,  high-turning,  first  stage  rotor  geometry.  It  is  scaled  up 
three  times  to  accommodate  the  cooling  scheme  and  instrumentation. 
The  span  is  15.3  cm  (6”)  and  the  aerodynamic  chord  is  13.6  cm  (5.4”). 
Pitch  and  axial  chord  are  1 1 .4  cm  (4.5”).  The  Reynolds  Number  based 
on  aerodynamic  chord  (13.6  cm)  and  unheated  exit  conditions  is  about 
6*  10^.  The  Mach  number  distribution  was  shown  to  correspond  to 
design  conditions  (Hale  et  al.,  1997). 
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To  simulate  the  shock  emerging  from  the  trailing  edge  of  a  NGV, 
a  shock  tube  creates  a  shock  wave  which  is  sent  through  flexible 
tubing  and  along  the  leading  edges  of  the  cascade  (see  Fig,  1).  The 
shock  strength  can  be  varied  to  obtain  realistic  pressure  ratios  by 
mounting  diaphragms  of  different  thickness  between  the  driver  and 
driven  section  of  the  shock  tube. 

Cooling  Design 

A  schematic  of  the  showerhead  film  cooling  design  is  shown  in 
Fig.  3  in  Part  I  (Popp  et  al.,  2000).  The  design  is  a  generic  showerhead 
layout  with  four  rows  of  holes  close  to  the  leading  edge  and  a  row  of 
“giir*  holes  on  both  pressure  and  suction  side.  The  rows  around  the 
leading  edge  are  normal  to  the  chordwise  tangent  and  inclined  30° 
with  respect  to  the  span  wise  direction.  The  “gill”  row  on  the  suction 
side  is  normal  to  the  spanwise  direction  and  inclined  30°  with  respect 
to  the  chordwise  tangent  while  the  “gill”  row  on  the  pressure  side 
emerges  at  an  angle  of  45°.  All  holes  are  fed  by  a  common  plenum 
which  is  insulated  to  reduce  internal  heat  transfer.  The  coolant  supply 
and  coolant  pressure  control  is  shown  in  Fig.  5  in  Part  I,  This  setup 
allows  the  control  of  the  ratio  of  coolant  total  pressure  to  freestream 
total  pressure.  This  corresponds  to  the  Momentum  Ratio,  since  air  is 
used  as  the  coolant.  The  objective  temperature  ratio  T/Tc  was  2,  but 
was  allowed  to  vary  during  the  run.  Therefore,  the  Density  Ratio  and 
Blowing  Ratio  also  vary  throughout  the  experiment,  while  the 
Momentum  Ratio  is  nearly  constant. 

Sensors 

The  six  measurement  locations  are  shown  in  Fig.  3  and  4  in  Part  I 
(Popp  et  al.,  2000).  At  each  measurement  location  a  triplet  of  gauges  is 
installed.  A  Kulite  Pressure  Sensor  XCQ-062-50a  is  monitoring  static 
surface  pressure.  A  Vatell  HFM-7/L  Heat  Flux  Microsensor  gives  a 
direct  reading  of  surface  heat  flux  and  of  surface  temperature,  A 
surface  thermocouple  is  used  as  a  calibration  tool  for  the  resistance 
temperature  device  of  the  HFM-7/L.  For  the  high  speed  measurements 
of  surface  pressure  and  heat  flux,  the  dynamic  characteristics  of  the 
heat  flux  and  pressure  sensors  had  to  be  known.  In  a  shock  tube 
experiment,  the  transfer  function  of  the  Kulite  XCQ-062-50a  pressure 
sensor  with  B -screen  was  found  to  be  modeled  well  as  a  second  order 
system  in  the  frequency  range  of  up  to  100  kHz.  The  transfer  function 
of  the  heat  flux  sensors  was  experimentally  determined  to  be  of  first 
order  with  a  time  constant  of  about  I7ps. 

Data  Acguisition  and  Signal  Processing 

Both  surface  static  pressure  and  heat  flux  were  recorded  at  500 
kHz  with  a  LeCroy  High  Speed  Waveform  Recorder.  The  signals  were 
filtered  at  102.4  kHz  using  Frequency  Devices  Model  9064  8-Pole 
Buiterworth  Low-Pass  Filters.  Since  the  transfer  functions  of  both 
pressure  and  heat  flux  sensors  show  significant  signal  attenuation  and 
phase  shift  in  the  frequency  range  below  the  cut-off  frequency,  the 
recorded  signals  were  corrected.  The  transfer  functions  that  had  been 
determined  experimentally  were  inversely  applied  to  the  pressure  and 
heat  flux  data  (Popp,  1999). 

To  find  the  mean  heat  flux  before  shock  impact,  a  time  window 
of  5  ps  before  the  sudden  rise  in  heat  flux  was  used  for  averaging. 

Test  Matrix 

A  technique  for  the  analysis  of  the  transient  heat  transfer  due  to 
shock  passing  was  presented  in  Popp  et  al.  (1999).  This  technique  will 
be  applied  to  heat  flux  time  histories  in  the  section  “Unsteady 
Decomposition”.  The  analysis  was  done  with  and  without  film  cooling 
and  under  a  wide  variety  of  parameters.  Without  film  cooling,  two 
parameters  were  varied:  Shock  strength  and  heat  flux  before  shock 


impact.  The  shock  strength  was  varied  in  two  steps.  In  terms  of  peak 
static  pressure  ratio  at  gauge  #2  the  recorded  peak  values  of  shock 
strength  were  1.2  and  1.3  when  the  pressure  signal  was  filtered  at  25 
kHz..  The  level  of  heat  flux  before  shock  impact  was  varied  by 
triggering  the  shock  at  different  times  during  the  tunnel  run.  As  seen  in 
Fig.  7  (a)  and  1 1  (a)  in  Part  I,  the  heat  flux  level  changes  during  the 
test.  It  was  attempted  to  obtain  three  different  levels  of  initial  heat  flux 
(high,  medium,  low).  The  two  values  of  shock  strength  and  the  three 
values  of  heat  flux  yield  a  test  matrix  of  six  experiments  without  film 
cooling.  For  the  tests  with  film  cooling  one  additional  paranwter  was 
varied.  The  ratio  of  cc^  ^ant  total  pressure  to  freestream  total  pressure 
was  set  to  the  nominal  value  (l.CW)  and  a  high  value  (1.2).  With  the 
shock  strength  and  t!ie  initial  level  of  heat  flux  varied  as  in  the  case 
without  film  cooling,  this  resulted  in  a  test  matrix  of  twelve 
experiments  with  film  cooling.  In  section  “Unsteady  Decomposition” 
only  sample  experiments  are  shown.  The  conclusions  drawn  there 
have  been  verified  for  all  conditions  tested  in  these  lest  matrices. 

In  sections  “Comparison  with  Rigby’s  Model”  and  “Comparison 
with  Moss’  Model”,  two  models  predicting  unsteady  heat  flux  from 
surface  static  pressure  are  compared  to  measured  heat  flux.  For  this 
purpose  pressure  and  heat  flux  time  histories  had  to  be  sampled  and 
filtered  at  the  same  frequencies,  and  the  signals  had  to  be  corrected  in 
order  to  have  the  same  signal  attenuation  and  phase.  This  was  done 
according  to  the  procedure  described  briefly  in  section  “Data 
Acquisition  and  Sighal  Processing”.  For  the  comparison  a  separate  and 
shorter  scries  of  experiments  was  done.  Pressure  and  heat  flux  were 
only  recorded  at  measurement  locations  #1  and  #2.  Two  tests  with  film 
cooling  and  one  test  without  film  cooling  were  done  at  the  higher 
shock  strength.  Since  the  frequency  content  of  the  pressure  data  is 
much  wider  in  these  tests  (102.4  kHz  compared  to  25  kHz),  the  peak 
unsteady  pressure  ratio  measured  at  gauge  #2  is  about  1.5  (compared 
to  1 .3  in  the  tests  with  the  lower  cut-off  frequency),  while  the  “real” 
shock  strength  is  identical.  For  the  two  experiments  with  film  cooling, 
the  ratio  of  coolant  total  pressure  to  freestream  total  pressure  was  set 
to  the  nominal  value  of  1.04.  One  test  with  and  one  without  film 
cooling  are  presented, 

RESULTS  AND  DISCUSSION 


General  Observations 

A  series  of  Shadowgraphs  taken  with  a  high-speed  camera 
showed  no  particularly  strong  interaction  between  the  shock  wave  and 


Figure  2:  Sample  Shadowgraph  Illustrating  the  Interaction  between 
Passing  Shock  and  Layer  of  Coolant  Air 
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the  layer  of  cooling  air.  No  film  detachment  or  separation  was 
observed.  Also  no  “Separation  Bubble”  nor  “Vortical  Bubble”  was 
seen.  Fig.  2  shows  a  sample  Shadowgraph  of  the  shock  passing 
process  with  the  shock  just  above  gauge  #1. 

Samples  of  the  heat  flux  traces  without  film  cooling  are  shown 
for  gauges  #1  through  #5  in  the  upper  graph  of  Fig.  3.  The  lower  graph 
in  the  same  figure  shows  the  pressure  traces  firom  gauges  #2,  #4,  #5 
and  #6  from  the  same  experiment.  Because  most  of  the  pressure  traces 
were  only  used  for  visual  comparison,  the  pressure  signals  shown  here 
were  filtered  at  only  25  kHz.  The  x-axes  on  the  two  graphs  are 
identical  to  illustrate  the  time  history  of  the  passing  shock  as  it  passes 
over  the  gauges.  The  y-axes  on  each  graph  are  also  identical  to  allow 
for  a  comparison  of  the  relative  strength  of  the  shock  on  each  gauge. 
The  time  mean  value  of  heat  transfer  before  shock  impact  was 
removed  to  show  only  the  unsteady  component  of  the  heat  flux  for 
better  comparison.  This  figure  illustrates  that  the  time  histories  of 
pressure  and  heat  flux  are  in  phase  and  similar  in  shape.  Higher  values 
of  unsteady  pressure  correspond  to  higher  values  of  unsteady  heat  flux. 
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Figure  3:  Traces  of  Unsteady  Heat  Flux  (Upper  Graph)  and  Unsteady 
Pressure  (Lower  Graph)  Due  to  Shock  Passing  Without  Film  Cooling 


The  same  observation  can  be  made  for  the  experiments  with  film 
cooling.  The  unsteady  pressure  and  heat  flux  histories  from  one 
experiment  with  film  cooling  are  shown  in  Fig,  4,  When  comparing 
Fig.  3  and  4  it  can  also  be  seen  that  the  traces  of  unsteady  pressure  and 
heat  flux  with  and  without  film  cooling  appear  similar  in  shape  and 


phase.  This  is  in  contradiction  to  the  statement  by  Rigby  et  al.  (1990), 
who  concluded  that  “the  unsteady  shock  wave  signature  changes  when 
film  cooling  is  introduced”.  For  the  front  suction  surface  they  observed 
similar  but  reduced  traces  of  unsteady  Nusselt  Number  with  the 
introduction  of  film  cooling.  This  observation  is  supported  by  Fig.  3 
and  4.  Even  though  the  pressure  traces  with  and  without  film  cooling 
are  nearly  identical  the  heat  flux  signatures  with  film  cooling  are  about 
20  to  30%  lower. 
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Figure  4:  Traces  of  Unsteady  Heat  Flux  (Upper  Graph)  and  Unsteady 
Pressure  (Lower  Graph)  Due  to  Shock  Passing  With  Film  Cooling 


“Unsteady  Decomposition” 

The  unsteady  component  of  heat  flux  due  to  shock  passing  can  be 
analyzed  by  a  general  decomposition  of  all  parameters  involved  (Popp 
et  al.  (1999)).  The  expression  for  the  unsteady  heat  flux  as  a  function 
of  all  time  dependent  parameters  is  then: 

q'  =  h'(T.,-Tj+hT:„+h'C  0) 


Die  superscript  '  refers  to  the  unsteady  component  of  the  physical 
parameter  while  the  variables  without  the  superscript  refer  to  the  value 
of  the  variable  before  shock  impact.  In  the  case  without  film  cooling 
the  adiabatic  wall  temperature  T^w  and  its  time  varying  component 
T",w  can  be  replac^  by  the  recovery  temperature  and  its 
corresponding  time  varying  ftinction. 


4 


Unsteady  Variation  of  Heat  Transfer  Coefficient.  The  first 
term  on  the  right  hand  side  of  Eq.  1  shows  that  the  contribution  of  the 
unsteady  heat  transfer  coefficient,  h',  to  the  overall  unsteady  heat  flux 
depends  on  the  temperature  difference  before  shock  impact,  or 

equivalendy  the  heat  transfer  before  shock  impact.  However,  it  will  be 
shown  in  this  section  that  the  overall  unsteady  heat  flux  is  independent 
of  the  value  of  heat  flux  before  shock  impact.  This  implies  that  the 
unsteady  component  of  heat  transfer  coefficient  is  not  contributing 
significantly  to  the  overall  unsteady  heat  flux.  If  h'  was  of  the  same 
order  of  magnitude  as  the  value  of  heat  transfer  coefficient  before 
shock  nnpact,  it  would  have  to  contribute  to  the  unsteady  heat  flux  by 
an  amount  of  the  order  of  the  value  of  heat  flux  before  shock  impact. 
This  is  clearly  not  the  case  as  seen  in  Fig.  5  and  6, 

Fig.  5  shows  the  superposition  of  heat  flux  data  from  three 
different  shock  passing  experiments  without  film  cooling  at  different 
levels  of  initial  heat  flux  (all  at  the  higher  shock  strength).  For  each  of 
the  five  gauges  the  upper  graph  shows  the  sum  of  the  steady  and 
unsteady  heat  flux  due  to  shock  passing.  The  different  levels  of  q 
before  shock  impact  can  be  found  in  these  graphs.  In  the  lower  graphs 
the  values  of  heat  flux  before  shock  impingement  were  removed, 
yielding  the  unsteady  component  of  heat  flux  q'.  Fig,  6  shows  the 
same  arrangement  for  three  experiments  with  film  cooling  (higher 
shock  strength,  coolant  to  freestream  total  pressure  ratio=1.2).  It  can 
be  seen  that  the  three  traces  of  unsteady  heat  flux  collapse  in  both 
cases.  This  illustrates  that  the  unsteady  component  of  beat  flux  is 
independent  of  the  initial  level  of  heat  flux  before  shock 
impingement,  thereby  proving  that  the  contribution  of  the 
unsteady  heat  transfer  coefficient  to  the  overall  unsteady  heat  flux 
is  small  both  with  and  without  film  cooling  (i.e.  h'  is  much  smaller 
than  h). 
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Figure  6:  Heat  Flux  and  Unsteady  Heat  Flux  Due  to  Shock  Passing. 
Upper  Graphs:  Overall  Unsteady  Heat  Flux  During  Shock  Impact. 
Lower  Graphs:  Unsteady  Component  of  Heat  Flux. 

Three  Experiments  with  Film  Cooling. 

(Shock  Pressure  Ratio  1.3,  Coolant  to  Freestream  Pressure  Ratio  1.2) 
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Figure  5:  Overall  and  Unsteady  Heat  Flux  Due  to  Shock  Passing. 
Upper  Graphs:  Overall  Unsteady  Heat  Flux  During  Shock  Impact. 
Lower  Graphs:  Unsteady  Component  of  Heat  Flux. 

Three  Experiments  without  Film  Cooling. 

(Shock  Pressure  Ratio  1.3) 


The  terms  containing  h'  in  Eq.  1  can  be  dropped  based  on  this 
conclusion.  For  the  experiments  without  film  cooling  this  yields  the 
following  relation  between  the  unsteady  variation  of  recovery 
temperature  and  unsteady  heat  flux: 


q'  =  hT; 


(2) 


For  the  experiments  with  film  cooling  the  corresponding  relation  is: 


q'  =  hcC 


or 


(3) 


These  equations  imply  that  all  the  unsteady  heat  flux  is  caused  by 
a  temperature  variation  induced  by  the  passing  shock,  which  follows 
from  the  conclusions  drawn  above. 

Unsteady  Variation  of  Film  Effectiveness.  The  effect  of 
the  passing  shock  on  the  film  effectiveness  can  be  analyzed  by 
decomposing  the  unsteady  variation  of  adiabatic  wall  temperature 
according  to  Popp  et  al.(1999): 

c  =t;  •  (1  (T,  (4) 


Or  rearranged: 

n'[(T,  -tJ+t;]=t;(i-ii)-c  (S) 
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The  difference  between  the  variation  of  the  recovery  temperature 
scaled  by  (1-t|)  and  the  variation  of  adiabatic  wall  temperature  is  a 
measure  for  the  magnitude  of  t|'.  The  scaling  factor  for  Ti'  on  the  left 
hand  side  of  Eq.  5  is  always  a  relatively  high  number  since  the 
difference  between  the  recovery  and  the  coolant  temperature  is  of  the 
order  of  150®C.  To  estimate  the  magnitude  of  “n',  the  two  terms  on  the 
right  hand  side  of  Eq.  5  are  plotted  in  Fig.  7.  T'r  is  calculated  from  an 
uncooled  experiment  shown  in  Fig.  5  using  Eq.2.  The  results  for  T'*. 
are  so  similar  for  the  three  experiments  that  a  separate  comparison  is 
not  necessary  and  a  representative  trace  of  T'r  can  be  used.  The  value 
of  q  used  for  each  gauge  location  is  die  local  mean  value  of  q  from  all 
three  experiments.  The  actual  values  were  determined  according  to  the 
procedure  explained  in  Part  I  of  this  paper  (Popp  et  al.,  2(XX)). 
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Figure  7:  Comparison  of  T'/l-q)  and  T'^^  for  Three  Different  Film 
Cooled  Experiments. 


It  can  be  seen  from  Fig.  7  that  the  difference  between  T^fl-q) 
and  T',w  does  not  exceed  7°C  at  all  times  in  all  three  experiments  and 
the  five  gauge  locations.  With  [(Tr-Tc)+T'J  on  the  order  of  150°C  this 
will  yield  values  of  q'  not  exceeding  5%.  This  value  is  relatively  small 
compared  to  the  value  of  q  (or.  the  order  of  20%  to  30%).  The 
variation  of  the  film  effectivercjs  can  therefore  be  considered  a 
secondary  effect  in  the  heat  transfer  process  of  shock  passage. 

Comparison  with  Riabv^s  Model 

In  1989  Rigby  et  al.  (1989)  published  a  modified  version  of  a 
model  for  the  prediction  of  unsteady  heat  transfer  first  developed  by 
Johnson  et  al.  (1988).  The  two  models  are  very  similar  but  the  later 
version  came  with  a  more  thorough  derivation.  Therefore,  Rigby’s 
model  will  be  used  for  comparison  with  the  data  presented  here.  The 
reader  is  referred  to  Rigby  et  al.  (1989)  for  the  detailed  derivation. 
Two  components  of  unsteady  heat  flux  are  identified.  One  of  them  is 
caused  by  compression  of  the  boundary  layer: 


qc=q- 


P  +  P 


(6) 


The  second  component  of  unsteady  heat  flux  is  caused  by 
transient  conduction  in  the  fluid  close  to  the  wall.  Neglecting 
convective  and  viscous  terms  in  the  boundary  layer  equation  it  is 
found  that  the  one-dimensional  heat  conduction  equation  applies  to 
this  component  of  surface  heat  flux.  The  temperature  time  history  used 
to  find  this  component  of  heat  flux  is  expressed  as: 


T(0,t)=T^ 


wV' 

p  J 


(7) 


This  temperature  time  history  is  processed  according  to  Oldfield  et  al 
(1978)  to  yield  the  second  component  of  unsteady  heat  flux  qj,; 

2  •  JkpCp  m-l  _ 

qh(t J= - i . +T„.,)-V;n-n  (8) 

V7C‘At  n=2 

The  overall  heat  flux  can  be  calculated  using: 

q  +  q'  =  qc+qh  (9) 


The  method  is  illustrated  in  Fig.  8  for  an  experiment  without  film 
cooling.  Only  data  from  gauges  #1  and  #2  are  available  for  this  type  of 
comparison.  The  graphs  on  the  left  hand  side  of  Fig.  8  show  data  from 
gauge  #1,  while  the  right  column  of  graphs  show  data  from  gauge  #2. 
The  top  graphs  show  the  pressure  time  history  and  the  time  history 
T(0,t)  calculated  according  to  Eq.  7.  The  second  row  of  graphs  show 
the  time  history  of  qh  and  qc  determined  according  to  Eqs.  8  and  6 
respectively.  The  last  row  of  graphs  shows  the  direct  comparison  of 
the  measured  and  predicted  heat  flux.  The  time  histories  of  measured 
heat  flux  and  surface  pressure  were  both  compensated  for  their 
specific  transfer  functions  in  the  frequency  range  below  102.4  kHz. 
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Figure  8:  Comparison  of  Measured  Heat  Flux  and  Heat  Flux  Predicted 
According  to  Rigby’s  Model  (Experiment  without  Film  Cooling). 
Left  Column:  Gauge  #1. 

Right  Column :Gauge  #2. 

Top  Row:  Time  Histories  of  Pressure  and  Th(0,t). 

Middle  Row:  qnandqc. 

Bottom  Row:  Comparison  of  Measured  and  Predicted  Heat  Rux. 
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From  the  second  row  of  graphs  it  is  clear  that  the  shock  effect  on 
qh  is  much  larger  than  on  qc  for  the  levels  of  q  and  the  pressure  ratios 
encountered  in  these  experiments.  The  last  row  of  graphs  showing  the 
direct  comparison  between  the  measured  and  predicted  heat  flux 
illustrates  that  the  model  overpredicts  the  unsteady  heat  flux  due  to 
shock  passing.  The  reason  for  this  must  be  that  at  least  one  of  the  three 
major  simpliHcations  made  in  the  derivation  of  the  model  was  not 
justified  (neglecting  convective  and  viscous  terms,  dropping  the  term 
containing  the  time  derivative  of  pressure  in  Eq.  (14)  of  Rigby  et  al. 
(1989)).  Fig.  9  shows  the  same  comparison  presented  in  Fig.  8  for  an 
experiment  with  film  cooling.  The  time  histories  of  pressure  are  nearly 
identical  with  and  without  film  cooling.  Therefore,  the  traces  of  Th(O.t) 
and  qh  are  also  similar,  and  qn  overpredicts  the  measured  heat  flux. 

Gauge  #1  Gauge  #2 


40 

P 

20'Z 


^10 

1. 


1 

—  (Rigby) 

—  (Rigby) 

L_ 

1.1  1.2 
Time  [me] 


1.3 


Time  [ms] 

Figure  9:  Comparison  of  Measured  Heat  Rux  and  Heat  Flux  Predicted 
According  to  Rigby’s  Model  (Experiment  with  Film  Cooling). 

Left  Column:  Gauge  #1 
Right  Column: Gauge  #2 

Top  Row:  Time  Histories  of  Pressure  and  Th(0,t) 

Middle  Row:  qh  and  q^ 

Bottom  Row:  Comparison  of  Measured  and  Predicted  Heat  Flux 


Comparison  with  Moss*  Model 

In  1995  Moss  ei  al.  (1995)  published  ensemble-averaged  heat 
flux  and  pressure  data  from  a  rotating  turbine  facility.  Using  Johnson’s 
model  -  which  is  almost  identical  to  Rigby’s  model  -  they  predicted 
heat  flux  from  the  pressure  lime  histories  and  found  reasonable 
agreement.  By  simply  assuming  a  constant  heat  transfer  coefficient 
and  predicting  the  heat  transfer  from  a  temperature  variation  induced 
by  isentropic  compression,  an  even  belter  prediction  could  be  found: 


qM=hTM  with  Tm=T- 


III 

^p  +  p'V 


-1 


(10) 


The  constant  temperature  T  in  Eq.  10  was  not  specified  by  Moss  et  al. 
(1995).  For  the  comparison  in  this  investigation,  the  freestream  static 


temperature  will  be  used.  Fig.  10  and  1 1  demonstrate  the  method  and 
compare  the  predicted  heat  flux  to  the  data  shown  in  the  last  section 
for  an  experiment  without  and  with  film  cooling  respectively. 

Gauge  «1  Gauge  «2 


Figure  10:  Comparison  of  Heat  Flux  Data  with  Moss’  Assumption  of  a 
Constant  h  (Experiment  Without  Rim  Cooling). 

Left  Column:  Gauge  #1. 

Right  Column:Gauge  #2. 

Top  Row:  Pressure  and  Temperature  Time  Histories 

Bottom  Row:  Comparison  of  Measured  and  Predicted  Heat  Rux 
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Figure  1 1 :  Comparison  of  Heat  Rux  Data  with  Moss’  Assumption  of  a 
Constant  Heat  Transfer  Coefficient  (Experiment  With  Film  (Tooling). 
Left  Column:  Gauge  #1 
Right  Column:  Gauge  #2 

Top  Row:  Pressure  and  Temperature  Time  Histories 

Bottom  Row:  Comparison  of  Predicted  and  Measured  Heat  Rux 
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Moss’  prediction  matches  the  measured  heat  flux  much  better  in 
terms  of  magnitude  and  shape  than  Rigby’s  model  presented  in  the  last 
section,  even  though  it  slightly  underpredicts  the  initial  peak.  The 
conclusion  from  section  “Unsteady  Decomposition”  is  supported  by 
this  observation.  The  contribution  of  the  variation  of  the  heat  transfer 
coefficient  to  the  overall  unsteady  heat  flux  proves  to  be  small 
compared  to  the  unsteady  heat  transfer  caused  by  the  variation  of 
temperature  induced  by  the  passing  shock. 

CONCLUSIONS 

The  unsteady  heat  transfer  caused  by  a  passing  shock  wave  on  the 
suction  side  of  a  transonic  turbine  blade  was  investigated  with  and 
without  film  cooling.  Heat  flux  data  was  analyzed  by  means  of  a 
mathematical  decomposition  and  by  comparison  with  predictive 
models  available  in  the  literature.  The  mathematical  and  comparative 
analysis  leads  to  the  conclusion  that  the  unsteady  heat  transfer  is 
mostly  driven  by  the  temperature  fluctuation  induced  by  the  passing 
shock  wave.  Heat  transfer  coefficient  and  film  effectiveness  are  not 
affected  enough  to  contribute  significantly  to  the  overall  unsteady  heat 
transfer.  This  conclusion  was  verified  on  six  locations  on  the  suction 
side  for  shocks  of  different  strength  and  for  different  coolant  blowing 
rates.  Since  the  temperature  in  an  actual  engine  is  fluctuating  around 
its  mean  value,  the  time  mean  heat  transfer  in  the  presence  of  shocks  is 
expected  not  to  differ  significantly  from  the  time  mean  heat  flux 
without  the  disturbance  of  the  shock. 

ACKNOWLEDGMENT 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific 
Research  (AFOSR)  under  grant  F0867 1-960 1062,  monitored  by  Dr, 
Jim  M.  McMichael,  Dr.  Mark  Glauser  and  Dr.  Tom  Beutner.  We 
would  like  to  thank  Messrs.  Scott  Hunter,  Monty  Shelton  and  Mark 
Pearson  of  General  Electric  Aircraft  Engines  for  their  collaboration  on 
this  project. 

REFERENCES 

Abhari,  R.S,  1996,”Impact  of  Rotor-Stator  Interaction  on  Turbine 
Blade  Film  Cooling”,  ASME  Journal  of  Turbomachinery,  Vol. 
118,  pp.  123-133 

Doorly,  D.J.,  Oldfield,  M.L.G.,  1985,*’Simulation  of  the  Effects  of 
Shock  Passing  on  a  Turbine  Rotor  Blade”,  ASME  85-GT-l  12 
Du,  H.,  Han,  J.-C.,  Ekkad,  S.V.,  1997,  “Effect  of  Unsteady  Wake  on 
Detailed  Heat  Transfer  Coefficient  and  Film  Effectiveness 
Distributions  for  a  Gas  Turbine  Blade”,  ASME  97-GT-166 
Dunn,  M.G.,  1989,  “Phase  and  Time-Resolved  Measurements  of 
Unsteady  Heat  Transfer  and  Pressure  in  a  Full-Stage  Rotating 
Turbine”,  ASME  89-GT-135 

Guenette,  G.R.,  Epstein,  A.H.,  Giles,  M.B.,  Haimes,  R.,  1989,  “Fully 
Scaled  Transonic  Turbine  Rotor  Heat  Transfer  Measurements”, 
ASME  Journal  of  Turbomachinery,  Vol.  1 1 1,  pp.  1-7 
Hale,  J.H.,  Diller,  T.E.,  Ng,  W.F.,  1997,  “Effects  of  Wake  on  Turbine 
Blade  Heat  Transfer  in  a  Transonic  Cascade”,  ASME  97-GT-130 
Heidmann,  J.  D.,  Lucci,  B.L.,  Reshotko,  E.,  1997,  “An  Experimental 
Study  of  the  Effect  of  Wake  Passing  on  Turbine  Blade  Film 
Cooling”,  ASME  97-GT-255 

Johnson,  A.B.,  Rigby,  MJ.,  Oldfield,  M.L,G.,  Ainsworth,  R.W., 
Oliver  M.J.,  1988,  “Surface  Heat  Transfer  Fluctuations  on  a 
Turbine  Rotor  Blade  Due  to  Upstream  Shock  Wave  Passing”, 
ASME  88-GT-172 

Johnson,  A.B.,  Oldfield,  M.L.G.,  Rigby,  MJ.,  Giles,  M.B.,  1990, 
“Nozzle  Guide  Vane  Shock  Wave  Propagation  and  Bifurcation  in 
a  Transonic  Turbine  Rotor’’,  ASME  90-GT-310 


Mehendale,  A.B.,  Han,  J.-C.,  Ou,  S.,  Lee,  C.P.,  1994,  “Unsteady  Wake 
Over  a  Linear  Turbine  Blade  Cascade  With  Air  and  C02  Film 
Injection:  Part  II  -  Effect  on  Film  Effectiveness  and  Heat 
Transfer  Distributions”,  ASME  Journal  of  Turbomachinery,  Vol. 
116,pp730-  737 

Moss,  R.W.,  Sheldrake,  C.D..  Ainsworth,  R.W.,  Smith,  A.D.,  Dancer, 
S.N.,  1995,  “Unsteady  Pressure  and  Heat  Transfer  Measurements 
on  a  Rotating  Blade  Surface  in  Transient  Flow  Facility”,  AGARD 
CP-571,  pp.  22-1-22-9 

Moss,  R.W.,  Ainsworth,  R.W.,  Garside  T.,  1997,  “Effects  of  Rotation 
on  Blade  Surface,  Heat  Transfer  An  Experimental  Investigation.”, 
ASME  97-GT-188 

Nix,  A.C.,  Reid,  T.,  Peabody,  H.,  Ng,  W.F.,  Diller,  T.E.,  Schetz,  J.A., 
1997,  “Effects  of  Shock  Wave  Passing  on  Turbine  Blade  Heat 
Transfer  in  a  Transonic  Cascade”,  AIAA-97-0160 

Oldfield,  M.L.G.,  Jones,  T.V.,  Schultz,  D.L.,  1978,  “On-Line 
Computer  for  Transient  Turbine  Cascade  Instrumentation,  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems,  Vol.  AES- 
14,  No,  5,  pp.  738-749 

Ou,  S.,  Han,  J.-C.,  Mehendale,  A.B.,  Lee,  C.P.,  1994, ’’Unsteady  Wake 
Over  a  Linear  Turbine  Blade  Cascade  With  Air  and  C02  Film 
Injection:  Part  I  —  Effect  on  Heat  Transfer  Coefficients”,  ASME 
Journal  of  Turbomachinery,  Vol.  1 16,  pp.  721-729 

Popp,  O.,  1999,  “Steady  and  Unsteady  Heat  Transfer  in  a  Film  Cooled 
Transonic  Turbine  Cascade,”  PhD  Dissertation,  Virginia  Tech. 

Popp,  O.,  Smith,  D.E.,  Bubb,  J.V.,  Grabowski,  H.C.  IB,  Diller,  T.E., 
Schetz,  J.A.,  Ng,  W.F.,  1999,  “Steady  and  Unsteady  Heat 
Transfer  in  a  Transonic  Film  Cooled  Turbine  Cascade”,  ASME 
99-GT-259 

Rigby,  M.J„  Johnson,  A.B.,  Oldfield,  M.L.G.,  Jones,  T.V., 
1989,’Temperature  Scaling  of  Turbine  Blade  Heat  Transfer  with 
and  without  Shock  Wave  Passing”,  International  Symposium  on 
Air  Breathing  Engines  ISABE  89-7070 

Rigby,  M.J.,  Johnson,  A.B.,  Oldfield,  M.L.G.,  1990,  “Gas  Turbine 
Rotor  Blade  Film  Cooling  With  and  Without  Simulated  NGV 
Shock  Waves  and  Wakes”,  ASME  90-GT-78 

Smith,  D.,  Bubb,  J.,  Popp,  O.,  Grabowski,  H.,  Diller,  T,  Schetz,  J., 
Ng,  W.,  2000,  “An  Investigation  of  Heat  Transfer  in  a  Film 
Cooled  Trar^  ^nic  Turbine  Cascade,  P^  I:  Steady  Heat 
Transfer”,  ASMa  2000-GT-202. 


8 


